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ABSTRACT 



We present models to predict high-frequency counts of extragalactic radio sources using physically grounded recipes to describe the 
complex spectral behaviour of blazars that dominate the mm-wave counts at bright flux densities. We show that simple power-law 
spectra are ruled out by high-frequency (v > 100 GHz) data. These data also strongly constrain models featuring the spectral breaks 
predicted by classical physical models for the synchrotron emission produced in jets of blazars (Blandford & Konigl 1979; Konigl 
1981). A model dealing with blazars as a single population is, at best, only marginally consistent with data coming from current 
surveys at high radio frequencies. Our most successful model assumes different distributions of break frequencies, vm, for BLLacs 
and flat-spectrum radio quasars (FSRQs). The former objects have substantially higher values of y^, implying that the synchrotron 
emission comes from more compact regions; therefore, a substantial increase of the BLLac fraction at high radio frequencies and at 
bright flux densities is predicted. Remarkably, our best model is able to give a very good fit to all the observed data on number counts 
and on distributions of spectral indices of extragalactic radio sources at frequencies above 5 and up to 220 GHz. Predictions for the 
forthcoming sub-mm blazar counts from Planck, at the highest HFI frequencies, and from Herschel surveys are also presented. 

Key words. Surveys: radio sources - - Radio sources: statistics - - Galaxies: active - - AGN: radio continuum 



1. Introduction 

In the last few years, large area surveys of extragalactic ra- 
dio sources (ERS) at high radio frequencies ^ 10 GHz - the 
Ryle telescope 15-GHz 9C s urvey {W aldram et all l2003h. the 
Very Small Array (VSA) one ("Watso n et al.L i2003l IClearv et al.L 
l2005h and, especially, the Australia Telescope AT20G surveys 
of the southern heni i sphere (ISadler et al.L l2006l iMassardi et al.L 
120081; iMurphv et al.', '2010") and the seven-year all-sky WMAP 
surveys (Gold et al., 2010) - have provided a great amount of 
new data on number counts, redshift distributions and emis- 
sion spectra of radio sources in the 10 - 100 GHz frequency 
range, which was poorly explored before. Moreover, number 
counts and related statistics of ERS are now available also at 
frequencies above 100 GHz, although they are estimated from 
smaller sky a reas and are co nsequently Um ited to fluxes 5 ^ IJy 
dVieira et al.L[2010; Maiiiage et al.Ll20TTb . 

At the beginning of the year 2011, available data on ERS 
have experienced an additional boost, by the publication of 
the Planck's Early Release Compact Source Catalogue (Planck 



ERCSC ) by the Planck Collaboration (Planck Collaboration 
201 lal) . In fact, the Planck ERCSC ( Planck Collaboration , 
201 Id) comprises highly reliable samples of hundreds of ERS at 
bright flux densities and high Galactic latitude, detected in many 
of the nine f requency channels (30-857 GHz) of the Planck 
satellite (Tau ber et al.L 12010). These samples of bright ERS - 
and the future ones that will be provided by the Planck Legacy 
Catalogue - will be unrivalled for many years to come at mm 
wavelengths. Indeed the statistically complete samples of ERS 
selected from the Planck ERCSC have already allowed the com- 
munity to achieve very relevant outcomes on number counts 
and on spectral properties in the 30-217 GHz frequency range 
(Planck Collaboration, 201 li). 

Early ev olutionary models of radio sources (e.g., 
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Da nese et all 1987; Dunlop & PeacockL Il990t iTofl^olatti et all 
19981: LTackson & Wall. 199^ " were able to give remarkable 
successful fits to the majority of data coming from surveys 
at y ^ 10 GHz, and down to flux densities of a few mJy. In 
particular, the model by Toffolatti et al. (1998) was capable to 
give a good fit of ERS number counts of WMAP sources as 
well, albeit with an offset of a factor of about 0.7 (see, e.g.. 
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iBennett et al.ll2003l) . Thus, this model was extensively exploited 
to estimate the radio source contamination of cosmic microwave 
background (CMB) maps ( Vielva et al., 2001, 2003; Tucci et al., 
12004 ) at cm/mm wavelengths. On the other hand, its very simple 
assumptions on the extrapolation of source spectra at mm 
wavelengths as well as new data published in the last ten years 
make it currently not up-of-date for more predictions, although 
it is still very useful for compari sons - even at v > 10 GHz - 
after a simple rescalin g (see, e.g.. Marriage et al.l 201 lb. 

Mo re recently, the De Zotti et al.l (l2005h and lMassardi et al.l 
(I2010ah new cosmological evolution models of radio sources 
have been able to give successful fits to the wealth of new 
available data on luminosity functions, multi-frequency source 
counts and also redshift distributions at frequencies ^ 5 GHz and 
^ 5 GHz, respectively. These two models are based on an accu- 
rate determination of the epoch-dependent luminosity functions 
(l.f.) for different source populations, usually separated by the 
value of the spectral index, a, of the observed emission spec- 
trum at low radio frequencies (1 ^ v S 5 GHz) by adopting a 
simple power-law approximation, i.e., S(v) oc v". This approxi- 
mation can be well assumed if there is synchrotron emission, but 
generally only in limited frequency intervals. 

If a > -0.5, sources are classified to have a "flat"-spectrum 
and are usually divided into flat-spectrum radio quasars (FSRQ) 
and BLLac objects, collectively called blazar^. Otherwise 
they are classified as "steep'-spectrum sources, which are 
mostly associat e d wit h powerful elliptical and SO galaxies 
(iToffolatti et all Il987h . Both populations consist of AGN- 
powered radio sources and the separation into two main source 
classes reflects that the region where the observed radio flux 
density is predominantly emitted is different in the two cases. 
For steep-spectrum sources, the flux originates in the extended 
(optically thin) radio lobes. For flat-spectrum sources the flux 
mainly comes from the compact (optically thick) regions of the 
radio jet. At lower flux densities, different populations of radio 
sources also contribute to the population of "steep'-spectrum 
radio sources: dwarf elliptical galaxi es, starburst galaxies , faint 
spiral and irregular galaxies, etc. (see iDe Zotti et al. I I2OO5I for a 
thorough discussion on the subject). 

The very recent, comprehensive review by De Zo tti et al.l 

(12010.) provides an up-to-date overview of all data published so 
far and at the same time of the cosmological evolution models 
and of the relevant emission processes that cause the observed 
ERS spectra. Moreover, an interesting and useful discussion of 
the relevant open questions on the subject is also presented. 

The predictions on high-frequency number counts of ERS 
provided by the above quoted cosmological evolution models 
are based on a statistical extrapolation of flux densities from the 
low-frequency data (< 5 GHz) at which the l.f.s are estimated. 
They adopt a simple power-law, characterized by an "average", 
fixed, spectral index, or by two spectral indices (at most) for each 
source population. Thus, this "classical" modelling has to be 
considered as a first - although successful - approximation, and 
it can give rise to an increasing mismatch with observed high- 
frequency (> 30 GHz) number counts. Indeed, the radio spectra 
in AGN cores can be quite different from a single power-law 



' Blazar sources are jet-dominated extragalactic objects - observed 
within a small angle of the jet axis - i n which the beamed com - 
ponent dominates the observed emission (I Angel & Stoc kmad [T980l) . 
They are characterized by a highly variable and polarized non-thermal 
synchrotron emission at GHz frequencies that is originated from elec- 
trons accelerated up to relativistic energies in the coll imated jets, which 
come out of the central active galactic nucleus (AGNi lUrrv & Padovanil 



in large frequency intervals. In particular, a clear steepening (or 
a bending down) at mm wavelengths is theoretically expected 
(Kellermann, 1966; Blandford & Konig], 1979). This steepen- 
ing has been alrea dy observed in some well known blazars 
dClegg et all 1 19831) and has also been statistically suggested by 
rece nt analyses of different ERS samples at v > 30 GH z (see, 
e.g., Waldram etal l l2007HGonzalez-Nuevo et"Zll2008h . 

The scenario giving rise to a moderate or to a more relevant 
(Aa > 0.5) spectral steepening in blazars - at frequencies of 
tens to hundreds GHz - is complex because different physical 
processes are intervening at the same time. The more relevant 
ones for our purposes are briefly sketched below (and will be 
discussed more extensively in Section 4). 

a) In the inner part of the AGN jet, i.e. typically at 
several thousand Schwarzschil d radii from the central col- 
lapsed object ( Ghisellini & Tave cchiol l2009h . various differ- 
ent and self-absorbed source components ("spherical blobs", 
very probably accelerate d by shock waves in the plasma; 
IBlandford & Koniglll 19791) are responsible of the observed syn- 
chrotron radiation at different distances from the central AGN 
core (Marscher& Gear, 1985; Marscher, 1996)0. In this situ- 
ation, the emerging jet- and blazar spectrum is approximately 
flat (a a 0.0, but it can be either very moderately steep, or even 
inverted); this is an indication that "the cores are partially opti- 
cally thick to synchrotron self- absorption, as expected in the jet 
model" (Marscher, 1996). 

b) However, the high-energy (relativistic) electrons 
- responsible for the self-absorbed synchrotron radi- 
ation, with the quoted a l most fl at emissio n spect rum 
jKeUermann & Paulinv-Toflil 119691: iMarscheit Il980al) - 
injected into the AGN jets are loosing their energy because of 
synchrotron losses, thus cooling down (i.e., electron ageing) 
when the rate of injection is not sufficient to balance the 
radiation losses. As a consequence, their synchrotron radiation 
spectrum has to steepen (at a given frequency) to a new a value, 
which depends on whether the e lectrons are continuously or 
instantaneously injected (see, e.g.. iKellermannL fl 966l Section II 
and Figure 5). 

c) In standard models of the synchrotron emission in blazar 
jets, the size of the optically thick core of the jet varies with 



freqtiency as rr oc v, ' ^* (see, e.g., KoiiiglL 19811 : IClegg et all 
119831: iLobanovi 120101: ISokolovskv et all 1201 Ih . "corresponding 
to the smallest radius, - rc, from which the optically thin 
synchrotron emission from the jet can be observed (see Sect. 
4)3 Therefore, the observed blazar sp ectrum at v > vm , now 
in the optically thin regime (see, e.g., IClegg et all Il983l) . will 
be moderately steep, with a typical a > ao + 0.5 value deter- 
mined by synchrotron-radiation losses owing to electron ageing 
( Kellermann. 1966 ) . 

As a result, at frequencies in the range 10 - 1000 GHz, de- 
pending on the parameters more relevant for the physical pro- 
cesses discussed here (see Sect. 4), the spectra of blazar so urces 
have to show a break or a turnover (see, e.g. Konigl, 119811 : 
IMarscheit \l99^ with a clear spectral steepening at higher fre- 
quencies. 



- Various examples of blazar spectra tentatively fitted by differ- 
ent synchrotron component s in the AGN jet are also given by 
iPlanck CollaborationI OOl lid) 

^ Recently, Very Long Baseline Array (VLB A) simultaneous obser- 
vations of the frequency-dependent shift seen in the core position of 
blazar sources (k nown as the "core shift") at nine frequencies in the 1.4- 
15.4 GHz range ( ISokolovskv et alil201lh have found new evidence of 
a parameter value k x I in the above quoted law, confiiTning theoretical 
predictions. 
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Moreover, in the frequency range where the CMB reaches 
its maximum, the interest in very accurate predictions of number 
counts of radio ERS is obviously increased, given that they con- 
stitute the most relevant contamina tion of CMB an i sotropy map s 
on small angular scales (see, e.g.. iToffolatti et all Il999l]2005h . 
Indeed, current high-resolution experiments at mm/submm 
wavelengths require not only to remove bright ERS from CMB 
anisotropy maps but also to precisely estimate the contribution 
of undetected point sources to CMB anisotropics. This problem 
is particularly relevant for current as well as forthcoming CMB 
polarization measurements (Tucci et al., 2004, 2005). 

In the present work we still compare our predictions on num- 
ber counts and other statistics of ERS to the most relevant ob- 
servational data sets at v > 5 GHz by extrapolating the spec- 
tral properties of the ERS observed at low radio frequencies. 
However, our approach is somewhat different than before. The 
fundamental difference consists in providing, for the first time, 
a statistical characterization of the sources' spectral behavior at 
high radio frequencies that takes into account - at least for the 
most relevant population of ERS at cm to mm wavelengths - the 
main physical mechanisms responsible of the emission. We fo- 
cus, in particular, on flat-spectrum radio sources, given that they 
are the dominant source p opulation in the flux- and freque ncy 
range we are interested in (iGiommi & Colafrancescoi l2004b . 

Finally, we analyse current data of another source popula- 
tion, although less relevant at bright flux densities: the popu- 
lation of the so-called "inverted'-spectrum sources (i.e., with 
a positive value of the spectral index a) detect ed in high- 
frequency radio surveys dballacasa et al., 2000; Bolton et al.L 
l200i ISadler etaLl 12006'; Murph v et al.ll2010l) . We give predic- 
tions also on their contribution to number counts and to other 
related statistics. In this latter case, our approach is still a purely 
statistical one, without dealing with the physical conditions that 
cause the observed emission. This simplified choice is justified 
by the current lack of data on the underlying physics of this 
source population and because these sources can be substantially 
contaminated by blazars obse rved in their active phase (see, e.g., 
iPlanck Collaborationi IMTlih . 

The outline of the paper is as follows: in Sect. 2 we dis- 
cuss number counts of flat- and steep-spectrum radio sources 
at 5 GHz; in Sect. 3 we analyse the spectral index distribu- 
tion of the different ERS populations; in Sect. 4 the basic as- 
sumptions of the simplified physical model for flat-spectrum 
sources is presented and extensively discussed; in Sect. 5 we sep- 
arately present the data on high-frequency peak spectrum (GPS) 
sources; in Sect. 6 we summarize the different model assump- 
tions for the various source populations we identify in Sects. 3, 
4, and 5; in Sect. 7 we present and discuss our model predictions 
on the extrapolation of the 5 GHz number counts and spectral 
properties of ERS to higher radio frequencies; finally, in Sect. 8 
we summarize our main conclusions. In addition, and for reduc- 
ing the main body of the article, we refer the reader to three ap- 
pendices: in Appendix A we give a brief description of the data 
sets currently available on ERS and analysed in this paper; the 
simplified formula adopted here for the estimation of the break 
frequency in blazars spectra is presented in Appendix B; finally, 
we discuss the main physical quantities that determine the value 
of the break frequency in Appendix C. 

2. Number counts of extragalactic radio sources at 
5 GHz 

The differential number counts of ERS at ~ 5 GHz are well 
known and have been extensively analysed and discussed (see 



iDe Zotti et all l2010l) . We plot in Fig.[T] the observed num- 
ber counts with the fits yielded by different model s, i.e., 
Toffolatti et al. ( 1998), De Zotti et al. (2005), and Massard ietal] 
dMOah . The plot covers the flux range between about 1 mJy 
and lOJy. At these levels the number counts are essentially 
dominated by AGNs, while at fainter fluxes the contribution of 
star-forming galaxies becomes increasingly important (see, e.g., 
Massardi et al., 2010a, and references therein for a thorough dis- 
cussion on the subject). 

However, the knowledge of the total counts of ERS at GHz 
frequencies is not sufficient for making predictions at higher ra- 
dio frequencies. Because of the presence of different source pop- 
ulations with different spectra, it is necessary to identify which 
populations dominate the number counts and their relative num- 
ber as a function of the flux density. To this aim, the large-area 
surveys of ERS available at GHz frequencies help us in classify- 
ing the observed ERS as steep- or flat-spectrum sourcesQ. 

In Fig.|2]we plot the differential number counts of these two 
source populations as calculated from different source samples, 
which are described in detail in AppendixlAl These results re- 
quire a careful discussion because the simple estimates of source 
spectral indices are not free from uncertainties, even in the case 
of negligible errors on published flux densities. First of all, we 
considered the catalogue of ERS with spectral information based 
on the NVSS and GB6 surveys (hereafter, the NVSS/GB6 sam- 
ple), which is defined in AppendixlAl Spectral indices from 
NVSS and GB6 data should be taken with some caution because 
of the different measurement epoch and the different resolution 
of the antennae. In particular, biased values of spectral indices 
could arise for resolution effects, especially when GB6 objects 
are resolved by the NVSS or have multiple components in the 
NVSS, which causes a flatter spectral index. On the other hand, 
variability mostly affects flat-spectrum sources. We discuss this 
specific problem in Section lTTI 

For estimating the magnitude of resolution effects on the 
4.8-GHz number counts of flat-spectrum sources, we identified 
all ERS resolved by the NVSS (i.e., with major axis > 45 arcsec) 
that are classified as flat-spectrum sources. We found 173 re- 
solved objects. Because flat-spectrum sources are usually as- 
sociated with compact objects with an emission dominated by 
the AGN core, these sources could be false identifications, and 
we marked them as steep-spectrum sources. Then we dealt with 
flat-spectrum sources with multiple NVSS counterparts and 
without clearly dominant NVSS objects: we looked for objects 
the brightest NVSS counterpart of which contributes < 75% 
of the total flux density and found 102. We redistributed these 
sources between the two populations according to the proportion 
of flat- and steep-spectrum sources in flux-density bins. After 
these corrections, the number counts of flat-spectrum sources 
were significantly reduced, especially at fluxes S ^ 300 mJy, 
whereas at higher fluxes the correction was small (see Fig.|2)|3 

Number counts of flat-spectrum sources in Fig. |2] were also 
obtained from the Parkes quarter-Jy sample (Jackson et al.i 
2002). The Parkes beam size is 8 arcmin at 2.7 GHz and 4 arcmin 
at 5 GHz, i.e., a factor 2 larger at the lower frequency. Therefore, 

Hereafter, we adopt the usual convention 5(v) oc y". Thus, we clas- 
sify an ERS as steep-spectrum if a < -0.5, and as flat-spectrum if 
a > -0.5, by using the flux densities measured at 1.4 and 4.8 GHz. 
Moreover, we consider as inverted-spectrum all those ERS for which 
a > 0.3. 

' These corrected counts will be not used later. However, they are an 
useful indication of the possible uncertainty in the spectral classification 
because of the different resolution between the two surveys discussed 
here. 
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Fig. 1. Differential number counts normalized to S^^^ for ERS 
at 4.8 GHz. Filled squares are observational data (see De Zotti 
et al. 2010). The plotted curves repr esent th e estimated number 
counts from the mod els by Toffolatti et al.l (11998,. dotted line), 
iDe Zotti et al] (I2005L blue dash-dotted line), and iMassardi et all 
(1201 Oai dashed line). Lower curves with the same line types as 
above represent the estimated number counts of flat-spectrum 
sources given by the same models. 



resolution effects should go in the opposite direction with re- 
spect to the NVSS/GB6 sample. Moreover, we increased the 
number counts of these sources by 10% (the percentage ob- 
served in the NVSS/GB6 sample) to take into account sources 
with -0.5 < a < -0.4, which were excluded from the original 
sample. 

Fig.|2l also show s the results from the CRATES catalogue 
(iHealev et all l2007h in the area of the GB6 survey. There, the 
number counts is computed taking into account 8.4-GHz mea- 
surements: we calculated the spectral index between 1.4 and 
8.4 GHz, and then the differential number counts at 4.8 GHz only 
for those sources that verified the condition o-j ^ > -0.5. The new 
source counts must be considered a lower limit because some re- 
ally flat-spectrum sources could have been discarded owing to 
variability or a spectral steepening between 4.8 and 8.4 GHz. We 
do not expect many of these, whereas the number of false iden- 
tifications coming from the classification of NVSS/GB6 data 
should be strongly reduced. If we compare these number counts 
with the ones derived from the NVSS/GB6 sample, we see that 
they partially agree at flux densities S ^ 0.5 Jy, whereas the dis- 
crepancy increases with lower fluxes. 

Finally, we also plotted the number counts of flat-spectrum 
sources estimated from the 5-GHz measurements present in the 
Austraha Telescope Compact Array 20 GHz survey (AT20G; 
iMumhv et al.[ l2010h : we limited our analysis to the almost- 
complete samples at declination 6 < -15° and flux lim- 
its 100 mJy and 50mJy, respectively (indicated as AT20G- 
dlSSlOO and AT20G-dl5S50, see AppendixEl; then, we took 
all sources with i) measurements at 5 and 8 GHz; ii) 5-GHz flux 
density S5 > lOOmJy; iii) > -0.5. The resulting number 



counts were corrected for the completeness of the samples and 
for the percentage of sources with 5- and 8-GHz data (89% 
and 84% respectively in the two AT20G sub-samples). We ob- 
serve that the number counts agree well with previous counts 
at S > 0.5 Jy, and with CRATES data at lower fluxes. ATCA 
measurements have the advantage to be nearly simultaneous and 
made with antennae of similar resolutions, which reduces there- 
fore the uncertainty in the spectral index estimates. Although not 
complete at 5 GHz, the sample built in this way provides reliable 
results starting from S <: 200 mJy, which should be considered 
at least as a lower limit for fluxes S < 1 Jy. 

In Fig.[TJ[2] the differential number counts predicted from 
the cosmolog i cal ev olutio n models of Toffo latt i et al. (19981) , 
iDe Zotti et all (l2005h . and lMassardi et al.1 (l2010ah are also dis- 
played. We see that whereas these models fit the total number 
counts at 5 GHz extremely well, their predictions on the num- 
ber of flat- or steep-spectrum ERS show significant differences 
and in general disag ree with the publis hed survey data. The 
iToffolatti et all ( 119981) and iDe Zotti et all ( l2005h models fit ob- 
servational counts of the two populations at flux S > 0.5 Jy quite 
well, but at lower fluxes they tend to overestimate (in the for- 
mer case) or to underestimate (in the latter case) the number of 
flat-spectrum sources. On the other han d, the ratio of steep- and 
flat-spectrum sources expected from the IMassardi et al.l(l2010ah 
model seems to be quite far from the observed ones, at least for 
5 > 0. 1 Jy. It is unlikely that uncertainties in the classification of 
radio sources could explain this great discrepancy. 

As a conclusion, at flux densities S > 100 mJy we can fit 
the differential number counts of flat-spectrum sources obtained 
from the observational data by a broken power law 



(1) 



where «(5 ) is the differential number counts normalized to S^^^. 
We find no = 47.4 Jy-'sr"', So = 1.67 Jy, = 0.50 and 
I - -0.66. We extrapolate this curve also to 5 < 100 mJy (rep- 
resented as a thick continuous line in the right panel of Fig.|2]l. 



3. Spectral indices of radio sources 

3. 1 . Spectral index distributions in tlie NVSS/GB6 sample 

It is usually assumed that the spectral index distribution of ERS 
at GHz frequencies can be described in all flux density ranges 
by the sum of two Gaussian distributions with maximum at 
a ~ -0.8 (for steep-spectru m sources) and at g ~ 0.0 (for flat- 
spectrum sources; see, e.g.. !De Zotti et al1l2010 V For checking 
this assumption, we calculated the spectral index distributions 
at 5 GHz of steep- and flat-sources in the NVSS/GB6 sample 
by dividing them into four flux-density intervals (see Table[T] 
and Fig.|3]l. On one hand, the distributions calculated for steep- 
spectrum sources seem not to change significantly with the flux 
density, peaking around -0.9 and having a median spectral index 
between -0.81 and -0.87. On the other hand, the distributions 
calculated for flat-spectrum sources show a clear dependence on 
the flux density interval, with steeper spectra at lower flux densi- 
ties. For S > 400 mJy the maximum in the distribution is around 
a ~ Q, which agrees w ell with observations from other surveys 
(e.g.. lRicci et alj|2004l) . However, at lower flux-density ranges, 
the maximum in the distributions gradually shifts to lower values 
of a and is around -0.4 for 5 < 251 mJy. The relative number 
of inverted-spectrum sources (a > 0.3) is about 14% at low flux 
densities, but it increases to almost 20% at 5 > 400 mJy. 




Fig. 2. Differential number counts at 4.8 GHz of steep-spectrum {left panel) and flat-spectrum {right panel) ERS from the 
NVSS/GB6 data (black empty squares). The predictions from the cosmological evolution models plotted in Fig. [1] are represented 
here by the same line - types as in tha t figure . In the left panel the thick continuous line represents the difference between the total 
number counts of the ' Tofi'olatti et al.l (119981) evolution model and our best fit to the observed number counts of flat-spectrum ERS 
(right panel; thick continuous line), i.e., it represents our best estimate of the differential number counts of steep-spectrum ERS. In 
the right panel we also plot the numb er counts of flat-sp ectrum ERS estimated from i) NVSS/GB6 data but corrected for resolution 
effects (black solid squares); ii) the J ackson et al.l (l2002h sample (red astericks); iii) the CRATES sub-sample in the GB6 sky area 
(blue empty triangles); iv) the 5-GHz sample of AT20G sources (cyan empty circles). See Sec|2]for more details. 

Table 1. Spectral properties of sources in the NVSS/GB6 sample. For each flux density range the columns indicate the total number 
of sources, the number of steep- and flat-spectrum (+inverted) sources, and the median values of their corresponding spectral 
index distributions; in the last column, the number of inverted-spectrum sources {a > 0.3) is also given. Cpeak and cr„ indicate the 
parameters of the best-fit truncated Gaussians for steep- and flat-spectrum sources (see Section 3.2). 
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3.2. Spectral index distributions vs. variability 

Because spectral indices are computed from non-simultaneous 
observations at 1.4 and 4.8 GHz, the source's variability could 
affect the spectral index distributions, especially for flat- 
spectrum sources that are observed to be variable at low frequen- 
cies as well. Accordingly, we expect that source variability will 
induce a higher dispersion in spectral index distributions and, re- 
lated to it, false identifications between steep- and flat-spectrum 
source populations. 

In iTingav et alJ (l2003h a study of source variability at GHz 
frequencies is carried out for 185 sources (167 flat- and 18 
steep-spectrum ERS): observations were made over a 3.5 year 
period with the ATCA at 1.4, 2.5, 4.8 and 8.6 GHz. For each 
source the authors provide the variability index, defined as the 
RMS fractional variation from the mean flux density. They found 
that at 4.8 GHz ~ 40% of flat-spectrum sources show a vari- 
ability > 10%, whereas only 5% vary by > 30%. Variability 
at 1.4 GHz is slightly lower On the other hand, no signifi- 
cant variability is observed in steep-spectrum sources, except 



for five sources (over 18) with a variability index between 4 
and 20 per cent. These results are consist ent with variability 
studies carried ou t at higher frequencies (ISadler et al.l l2006t 
iBolton et al.Ll2006h . 

We assume that flux densities vary at 1 .4 and 5 GHz ac- 
cording to the Tingay et al. results, which we summarize in 
Table|2] We generated a simulated sample of steep- and flat- 
spectrum sources at 5 GHz with a flux density distribution con- 
sistent with the GB6 number counts and with spectral index dis- 
tributions given by Gaussian functions truncated at ff = -0.5 
(with peak position o-peak and dispersion cr„). The 1.4-GHz flux 
density is calculated as 55(1.4/5)". After introducing variabil- 
ity in flux densities at 1 .4 and 5 GHz, we computed the new a- 
distributions by taking into account the uncertainties in GB6 and 
NVSS flux densities as well. Finally, we found the parameters 
of the truncated Gaussian distributions (o-peak and cr„), for which 
the simulated ff-distributions give the best fits of the observed 
NVSS/GB6 distributions. The results are reported in TablejTland 
in Fig. [3] where we show both the original truncated Gaussians 
and the final a-distributions produced from them. The value of 
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Fig. 3. Distributions of the spectral index calculated from 
the NVSS/GB6 sample for different flux density inter- 
vals: /i) S =[100,158]mJy; 12) S =[158,251] mJy; 13) 
S =[251,400] mJy; 14) S >400mJy. The thick black lines are 
the "best-fit" truncated Gaussians, and the thin red lines are the 
total distributions obtained from them after introducing source 
variability (see the text). 

Table 2. VariabiHty indices at 1.4 and 4.8 GHz for flat- and 
steep-spectrum sources. Columns 2, 3, and 5 give the frac- 
tions of sources that show the corresponding variability index 
in columns 1 and 4. For steep-spectrum sources there is no sig- 
nificant difference between the variability indices at 1 .4 and 4.8 
GHz. The estimated values are from iTingav et al.l (l2003l see 
text). 
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0.23 


4-20% 


0.30 


20-30% 


0.07 
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Q'peak is -0. 14 for S > 400 mJy, and steadily decreases to -0.4 in 
the lowest flux-density range. We see that simulations can repro- 
duce the NVSS/GB6 results very well. Variability only partially 
affects the spectral index distributions, increasing their disper- 
sion, and makes some flat-spectrum sources move to the steep- 
spectrum class. This effect is particularly relevant for sources 
with spectral index close to -0.5 and at flux densities < 250 mJy. 

The previous simulations also give us the number of sources 
that are misclassified because of variability. We found that the 
number of flat-spectrum sources from the NVSS/GB6 sample 
can be underestimated by about 5-8 per cent, whereas the num- 
ber of inverted-spectrum sources can be overestimated by only a 
few percent. Owing to their very low variability, the number of 
steep-spectrum ERS classified as flat-spectrum sources is small, 
and only down to S ~ lOOmJy can partially compensate errors 
in the classification of flat-spectrum ERS. 
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Fig. 4. Left panel: Distributions of 1-5 GHz (solid line) and 5- 
20 GHz (dashed line) spectral indices of steep-spectrum sources 
selected from the AT20G catalogue. Right panel: distribution of 
the difference a? - a?. The thin line is the Gaussian fit. 



3.3. High-frequency Steepening of steep-spectrum sources 

A high-frequency spectral steepening is expected in steep- 
spectrum sources due to electron ageing (iKeUermannl Il966h 
and has been observed in multifrequency surveys a t v > 5 GHz 
(iBolton et al.L 2004; Ric ci et alil2006l) . lRicci etaP ((2006) com- 
pared the 2.7-5 GHz and 5-18.5 GHz spectral indices and found 
a median steepening of Aa - 0.32. 

A spectral steepening in this frequency range is also ob- 
served in the AT20G data: we selected sources in AT20G- 
dl5S50 with 5-GHz flux density S > 500 mJy and spectral 
index ff j < -0.5. We found 211 objects. In Fig.|4]we plot the 
spectral index distribution computed between 1-5 GHz and 5- 
20 GHz, and the difference Aa = o-j - of . A median steepen ing 
of 0.28 is found, in agreement with the Ricci et al ] (12006') re- 
sult. The distribution of Aa is well described by a Gaussian with 
< Aa >= 0.28 and cr ^ 0.20, except for the negative tail. Values 



of Aa < (see also the presence of sources with af > -0.5 in 
the left pan el of Fig.lHl can arise both from sources with spectra 
that flatten dTucci et al.L 12008 ) and also from source variability 
because 1 GHz data are not simultaneous with the AT20G mea- 
surements. 

At frequencies v > 20 GHz a further spectral steepening is 
still possible: a recent indicati on in this sense is com ing from the 
follow-up of 9C sources (see lWaldram et aLll2007h . 



3.4. Higii-frequency spectral beiiaviour of flat-spectrum 
sources 

We focus now on the spectral behaviour of flat-spectrum sources 
at frequencies v > 5 GHz. The AT20G survey, thanks to its 
near-simultaneous observations at 5, 8 and 20 GHz, allows 
us to study in a quite broad range of frequencies spectral 
shapes that sho uld be not significantly affected by variability 
dMassardi et al.L BoiObi) . To this aim, we selected the sources 
of the AT20G-dl5S50 sample with S5 > 200 mJy and with 
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5-8 GHz spectral index > -0.5. This sample of 798 ob- 
jects is about 90% complete for sources with o-g" ^ -0.5, but 
it can lose objects that undergo spectral steepening between 8 
and 20 GHz. In Fig. |5] we show the distributions of and a^^ 
for those objects. The distribution of a^^ is shifted to lower val- 
ues by about 0.1/0.2 with respect to the distribution, with a 
large tail of objects with steep spectral index, extending up to 
-1.5. The percentage of sources with o-?" < -0.5 is ~ 22%. The 



median of distributions, excluding inverted-spectrum (a^ > 0.3) 
sources, is -0.08 and -0.27 for 5-8 GHz and 8-20 GHz spec- 
tral indices, respectively. Fig. |5] also shows the histogram of 



In this case we separate sources with 5-8 GHz 
spectral indices that are moderately steep (-0.5 < < -0.1) 
and flat (-0.1 < Cj < 0.3). We found that sources with flatter 
5-8 GHz spectral index tend to have a more relevant steepening 
between 8 and 20 GHz (the medians Aa are 0.13 and 0.23, re- 
spectively). An interesting feature in both distributions is also 
the presence of sources with typical steepening around 0.7-0.9, 
which seem to produce a secondary peak in the Aa distribution 
of sources with -0.5 < < -0.1. Because of the incomplete- 
ness of the sample for objects with ofg" < -0.5, the peak could 
be more prominent than observed, and it could extend to Aa ^ 1 . 
Finally, we note that most inverted-spectrum sources between 5 
and 8 GHz (o-^ > 0.3) have af < 0.3, and only ~ 12% ai-e still 

inverted up to 20 GHz. 

Our findings agree well with the results of iMassardi et al.l 
(1201 01^). A very strong average spectral steepening has also been 
found in the brightPACO sample (Massardi et al., 2010c), where 
the typical differences between the 5-10 GHz and 30^0 GHz 
spectral indices are in the range 0.5-1. 

Additional indications of a spectral steepening in flat- 
spectrum sources at frequencies higher than 20 GHz are given by 
OCRA-p observa tions at 3 1 GHz of CRATES selected sources 
(iPeel et al.Ll20Tqh and bv ATCA ob servations at 95 GHz of 130 

TTsi 



AT20G sources (Sadler et al.', '2008*). These results are also sup- 
ported by WMAP data (Gonzalez-Nuevo et al., 2008]at Jy lev- 
els, and by ACT data at 148 GHz dMarriage et aUfeOlTI) 



4. A spectral model for flat-spectrum sources 

As discussed in the previous section, recent multi-frequency sur- 
veys have clearly shown that ERS with flat spectra at GHz fre- 
quencies can present a downwards spectral curvature already at 
few tens of GHz. This spectral behaviour could well be an in- 
dication that the transition in the synchrotron spectra from the 
optically-thick to the optically-thin regime is observed, and that 
this transition can occur even at cm wavelengths. This conclu- 
sion is also supported by the results on the compactness of radio 
sources as a function of the spectral index, obtained from the 
AT20G survey. As an example, Massardi et al. (2010b, see their 
Fig. 4) do find a remarkable separation between steep-spectrum 
(q-j < -0.5) "extended" sources and flat-spectrum (o-j > -0.5) 
"compact" sources. However, if the compactness is considered 
versus the spectral index between 8 and 20 GHz, there appears to 
be also a population of compact and steep-spectrum ERS, corre- 
sponding to the 12% of the sample, which is not present at lower 
frequencies and suggests a spectral steepening at v <20GHz. 

By following this hypothesis, the key point for making pre- 
dictions on the high-frequency contribution of flat-spectrum 
ERS to number counts is to make an estimate of the frequency 
at which the spectrum of this source population can experience a 
substantial steepening (Aa ^ 0.5). To this aim, we use a simple 
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Fig. 5. Upper panel: distributions of spectral indices a^ (solid 
line) and (dashed line) for sources selected in AT20G- 
dl5S50 with ^5 > 200mJy and > -0.5. The dotted fine 

for inverted-spectrum sources 



shows the distribution of 



(a^ > 0.3). Lower panel: distribution of Aa for sources with 



line). 



-0.1 (solid Une) and with 



-0.1 <al< 0.3 (dashed 



physical description of the AGN radiation emission in the inner 
part of their relativistic jets. Our predictions are then calculated 
on the basis of typical physical parameters in current models for 
the emission of blazar sources. 
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4.1. Break frequency in synchrotron spectra from relativistic 
blazars jets 

It is generally believed that flat spectra observed in radio 
sources result from the supeiposition of different compo- 
nents of relativistic jets, each with a different synchrotron 
self-absorption fr equency (IKellermann & Pauliny-TodH 1 19691: 
ICotton et al.L[T98Q !). Relativistic shocks in the jet produce flares 
that are responsible for radio flux density variations on in- 
tervals from month s to years (e.g., Marscher & Gear, 19851 
IValtaoia et aLlll992l) . 

Flat radio spectra are predicted by models of synchrotron 
emission from inhomogeneous unresol ved relativistic jets 
(iBlandford & Konig]Lll979l:lKonigl . rT98Th . According to them, 
the observed flux density at a given frequency is dominated 
by the synchrotron-emitting component, which becomes self- 
absorbed at that frequency, because synchrotron self-absorption 
is the most likely absorption mechanism acting in the compact 
conical jet ( Konigli 119811) . When moving to high frequencies, 
self-absorbed synchrotron emissions from more and more com- 
pact regions are thus observed closer and closer to the AGNs 
core. However, these models also predict that at a particular fre- 
quency, jet emissions cease to be dominated by optically-thick 
synchrotron emission because of electrons cooling effects, and 
become dominated by synchrotron emission in optically-thin 
regime. Spectra then change from "fla t" to "st e ep". 

Let us review in more detail the iKonigll (Il981h model. In 
this model, the magnetic field H and the distribution of rela- 
tivistic electrons (described by a power law n(y) = /Ty"''"^"*, 
where y is the Lorentz factor) are considered inhomogeneous 
and scale with the jet radius as H{r) oc r^'" and K(r) oc r^". We 
assume m = 1 and n = 2, which coiTesponds to assuming the 
magnetic and the electron energy densities are in equipartition 
(IBlandford & Konigl. 1979). 

For the wavelengths range we are interested in, the syn- 
chrotron spectra from a portion of a relativistic conical jet at a 
distance r from the AGN core show two characteristic frequen- 



- v.5,„, the frequency where the synchrotron spectrum has the 
maximum and is estimated by setting the optical depth equal 
to unity. Below Vj„, the emitting source is optically thick to 
synchrotron self-absorption and the flux density rises as v^^, 
whereas at v > v„„ the emission has the optically-thin syn- 
chrotron spectrum 5 oc v" witho' < -0.5. Under the previous 
hypothesis, the maximum frequency scales as v,,,, oc r"'. It 
means that the size of the optically thick "core" of the jet 
decreases as the observed frequency increases. 

- Vsb (> Vsm), above which spectra steepen owing t o 
synchrotron-radiation losses (IBlandford & Konigli 1 19791) . 
The cutoff occurs when the synchrotron cooling time equals 
the reacceleration time of electrons. Contrary to v,,,,, Vsh in- 
creases as the distance from the core increases, proportion- 
ally to r. This implies that in an AGN jet there is a radius 
rM at which Vs„,{rM) - Vsh(rM). Thus, is the smallest 
radius from which optically-thin synchrotron emission can 
be observed with its characteristic "steep" spectral index, 
a < -0.5. 

As said before, the flux density observed from an unresolved 
conical jet is given by the sum of the emission from the different 
portions of the jet. The flat part of synchrotron spectra is ob- 
served up to the frequency vm (that we call "break frequency"), 
i.e. the frequency at which Vsm(rM) = VstirM)- At frequencies 
V > Vm, contributions from regions r < km inside the jet are 



negligible owing to the cooling of high-energy electrons. The 
observed flux density is thus dominated by optically-thin emis- 
sion from radii r > (v/vM)rM that verify v^hir) > v (Konigl 

Therefore, spectra at cm/mm wavelengths can be approxi- 
mated by two power laws 



5(v) = 



S{VM)iv/VMrf' 

S(vm)(v/vm)"" 



if V < Vm 
if V > Vm , 



(2) 



where aji and a,, represent the effect of the nonuniformity of 
sources in the optically thick and optically thin regimes. 

The frequency vm an d t he radius rM are provided by 
Blandford & Konigl (1979) and Konigl (1981) as the function of 
the relevant physical quantities of AGN jets. These expressions 
are quite complex, however, and depend on too many model pa- 
rameters. For simp licity (and, e .g., in agreement with iMarscheil 
1987 and GhisellinieLaljliaaJ, we assume that the flux density 
observed at frequency ^ vm is dominated only by the emission 
from the region at radius rM, and that other contributions are 
negligible. Then, for single emitting regions, synchrotron spectra 
can be approximately described by the homogeneous spherical 
model, and the break frequency is given by (Pacholczvk, 1970) 



Vm S ^ ' H ' (1 + z) ' 6 ' 



(3) 



where is the observed flux density at vm and 6 is the Doppler 
factor. Finally, 6 is the observed angular dimension of the emit- 
ting region: for a narrow conical jet of semiangle 0, the axis of 
which makes an angle ip with the dkection of the observer, 6 can 
be written as 



0: 



(l+zf (i+zf 
2 — rM tan(0) cos i// ^ 2 — ;:; ( 



D, 



(4) 



where Di is the luminosity distance of the source, and we as- 
sume the view a ngle i/r ^ 10° and cosiA ^ 1. As a typical valu e, 
we take 4> ~ 0.1 ( Konigl. 1 98 iVGhisel lini & Taveccliioll2009h . 

Under the hypothesis of equipartition condition, we can also 
have an estimate of the intensity of the magnetic field in the 
emitting region. Indeed, from the minimum energy conditions 
(Pacholczyk, 1970 ), the magnetic field intensity is 



24 



V 7 

where the minimum energy density Mmin is 



.10-1^ 



4/7 



(5) 



(6) 



and L is the luminosity of the emitting component, V the vol- 
ume in pc^. In Eq. ^ we are assuming that the total energy of 
electrons is equal to the total energy of protons. 

In the literature there is no clear agreement on whether the 
condition of equipartition is valid in blazar jetfl The equiparti- 
tion magnetic field should be considered only as an upper limit. 



* According to iMaraschi & Tavecchid (l2003h . the condition is veri- 
fied in flat-spectrum quasars, but not in BL Lacs, where the magnetic 
field energy is observed to be lower than the relativistic particle energy. 
Indications that the brightness temperature of the radio emission in jet 
components rapidly drops to the equipartition limit moving away from 
the core are also given by Lobanov (2010) and references therein. On 
the other hand, Ghise llinj et al . (2010) argue that blazar jets should be 
matter-dominated on scales where most of their luminosity is produced, 
and sub-equipartition magnetic fields are found only on smaller scales. 
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Fig. 6. Distribution of the break frequency, vm, for 10^ sources 
with 5-GHz flux density of 1 Jy (soHd Une), 0. 1 Jy (dashed Une) 
and lOmJy (dotted Hne). In this plot we only consider 0.01 < 
rM < lOpc (case CI). 

but because of the weak dependence of vm on H (vm <^ H^^^), 
we can safely adopt the equipartition approximation. 

The total luminosity in a homogeneous spherical source is 
related to the break frequency and to the flux density at vm as 
(the proof is given in Appendix B) 



Dl 6"-^ 



(l+z) 



(7) 



The factor /(a) depends on the optically-thin spectral index and 
on the lower and upper cutoff in the electron distribution func- 
tion. 

Finally, if the flux density at 5 GHz is known, Sm can be 
extrapolated according to Eq. ([2), S m - SsivM/Sy^'. By substi- 
tuting in Eq. Q the equipartition magnetic field //eq calculated 
above, the break frequency vm then only depends on the phys- 
ical parameters z, 6 and rM, and on the spectral indices before 
and after the break frequency, vm- Equation (O thus becomes 



C(a,a/,)[///(l+zy^-5/'">t]'^^ 



(8) 



withyS ~ 1 +0.5afi,/3D ~ 1,/?, ~ -1.5, /3s ~ -0.5 and/?, ^ -1. 
The explicit derivation of this equation is given in Appendix B. 

4.2. Break frequency and magnetic fields in fiat-spectrum 
sources 

The most relevant parameter for the calculation of the break fre- 
quency in blazars spectra is, undoubtedly, the radius r^. 

The radius is the distance from the AGN core of the jet 
portion that dominates the emission at vm, and it is related to 
the break frequency approximately as vm <^ (in agreement 
with predictions of inhomogeneous models). This parameter is 
the most relevant one in the estimate of vm because it defines the 
dimension and, thus, the compactness of the emitting region at 
that frequency. At the same time it is also the most critical one, 
because of the current uncertainty on its real value and of the 
relatively great range of possible values. 
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Fig. 7. Distribution of the break frequency, vm, for lO'* sources 
with 5-GHz flux density of 10/1/0. 1/O.OlJy (see labels in the 
Figure); the values correspond to the cases C2Co (black 
lines) and C2Ex (blue lines). In the two upper panels only model 
predictions corresponding to the flux densities of 0.01 and 0. 1 Jy 
are shown; in the two lower panels, only predictions correspond- 
ing to 1 and 10 Jy are shown. In the lower panel on the left side 
we also plo t the vm (red vertical arrow) of 3C 273 estimated by 
IClegg etal.l(ll983h . 



In the radio, i.e. synchrotron, regime the AGN jets emis- 
sion is expected to be produced at distances along the jet 
starting from 10~^-10~' pc up to parsec s cales (Marschei 
1996| ILobanovl 119981; iMaraschi & Tavecchid [20031; iLobanov , 
201(|). Moreover, the standard models used to fit the spec- 
tral energy distribution (SED) of blazars - i.e. the leptonic 



the l e ptonic 

ho mogeneous one-zone mo del (e.g. iGhiselUni et al.L 119981; 
Ghi selUni & Tavecchid l2009l) - typically find that the "dissipa- 
tion region" is at a distance from the black hole of around 10"^- 
10 'pc (Ghisellini &Tavecchio, 2009; Ghisellini et al., 20ld). 
However, the parameters of one-zone models are chosen to fit 
the high-energy emission of blazars (v ^ lO'"' Hz), whereas they 
are not able to reproduce the blazar spectra at cm/mm wave- 
lengths, where the emission is supposed to come from more ex- 
ternal and, thus, larger regions of the jet. Therefore, we consider 
the above quoted values as lower limits for r^. Emission from 
ultra-compact jet components are also required for explaining 
obser ved spectra that keep flat up to frequencies of ^ 10'^ Hz 
(e.g., lAbdo et al.L 120101; iGonzale z-Nuevo et all 12010*). On the 
other hand, observed breaks in blazars spectra at few tens of 
GHz need higher values of r^?. Upper limits for rM are pro- 
vided by VLBI observations, the resolution of which is of frac- 
t ions of milliarcseconds at frequencies v S 5 GHz (see, e.g., 
Sokolovs kv et all 1201 1 ^. The typical linear dimensions of the 
dominant AGN jet components derive d from these observa- 
tions are of the order of a few pa rsecs (iGhiselUni et al.l 119931 ; 
iJiang et all [19981 ILobanovl I20T0I) . 

To summarize, the range of possible values for should be 
0.01 < rM < 10 pc (i.e., 3 x 10'^ < tm < 3 x 10''^cm). This is 
a very large interval and it is taken as a first working case to cal- 
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culate the break frequency (hereafter, case CI). In addition, we 
consider more restricted intervals of km values, distinguishing 
between FSRQs and BL Lac objects: we take 0.01 < tm < 0.3 pc 
for BL Lacs and two different cases for FSRQs, 0.03 < < 1 pc 
(hereafter, C2Co, i.e. "more compact") and 0.3 < tm < 10 pc, 
a factor of 10 larger (hereafter C2Ex, i.e. "more extended"), 
values are assumed log-uniformly distributed inside the above 
ranges. This separation in these two classes of blazars is based 
on the fact that they show different spectral energy distribu- 
tions SEDs. On one hand, FSRQs, which are more powerful ob- 
jects, are observed to have a synchrotron peak frequency, VpW at 
lower frequencies, around lO'^-lO''* Hz, without a clear corre- 
lation with the radio power On the other hand, in BL Lacs Vp is 
found to be at higher frequencies and with a larger range of val- 
ues, from 10'^^ up to > 10'^ Hz (Foss ati et al.lll998tlAbdo et all 
I2OIO). This difference is motivated by the fact that BL Lacs are 
characterized by a lower intrinsic power and by a weaker exter- 
nal radiation field. Consequently, in BL Lacs the cooling through 
radiation losses is less dramatic and the electrons can be present 
with sufficiently high energies to still mantain synchrotron emis- 
sion up to these high frequencies (Ghisellini et al., 1998). That 
is why we assume that BL Lacs are in general more compact ob- 
jects compared to FSRQs, and that the emitting region at the 
break frequency, vm, is closer to the AGN core. As a conse- 
quence, the break frequency will be found to be on average at 
higher frequencies in BL Lacs than in FSRQs. 

It is important, however, to stress that there is no physical 
clearly established relationship between the break frequency vm 
and the peak frequency Vp in the synchrotron SED of AGNs. 
Their relative position can vary a lot from source to source. 
Indeed, although blazar spectra above vm are dominated by 
synchrotron emission in the optically-thin regime, the SED of 
blazars can still continue to increase up to frequencies » vm 
if the spectral slope is -1 ^ a ^ -0.5. Interesting exam- 
ples of SEDs for a complete sample of 104 northern blazars 
with average fluxes S > 1 J y at 37 GHz are presented by the 
iPlanck Collaboration! (|201 Ikl) . from which it is possible to have 
an idea of the very different positions of vm and Vp in the ob- 
served SEDs. 

Apart from r^, the other physical quantities relevant for cal- 
culating Vm from Eq. dS) are the redshift distribution, the Doppler 
factor 6 in AGNs jet, and the spectral indices before and after the 
break frequency, vm- Values of these physical quantities are dis- 
cussed in AppendixICland in Section|6l 

Figures|6]-|7] show the distribution of break frequencies in 
the observer frame for sources with the 5-GHz flux density of 
55 = 10, 1, 0.1 and 0.01 Jy, obtained by Eq.[8]for the different 
choices of discussed above. In general, we see that vm is on 
average lower for fainter sources, as expected from the relation 

°^ S'^M Eq. (O. for CI the large interval of possible val- 
ues of rM is reflected in an almost flat distribution of the break 
frequencies in the range 20-lOOOGHz, 10-600GHz, and 2- 
200 GHz for 55 = 1, 0.1 and 0.01 Jy, respectively. A few percent 
of sources with break frequency values vm ^ 10'^ Hz are also 
found. On the other hand, when conditions on are stricter, 
the distributions of vm are narrower In the cases C2Co-C2Ex 
we find that most of BLLacs of 10-1 Jy have a spectral break 
at V » 100 GHz, with a significant fraction at v ^ 10^ GHz. 
Fainter BL Lacs (100-lOmJy) have v^-distributions peaking 
around 50-200 GHz. Therefore, according to our results, BL 
Lacs should be observed at cm/mm wavelengths with almost flat 



' Vp represents the frequency at which a maximum is reached in the 
synchrotron SED of AGNs (in terms of vF,,). 




H., (mG) 

Fig. 8. Distribution of the magnetic field strength in equipartition 
conditions at a distance of 1 pc from the AGN core. Black curves 
are for FSRQs and red curves for BLLacs; solid curves are for 
sources with 5 5 — 1 Jy and dashed curves if ^5 = 0.1 Jy. 



spectra, at least for flux densities > 0.1 Jy. The vm distributions 
of Jy BL Lacs extend up to 10'"^ Hz, as does the one for CI. 

For FSRQs, the break frequency varies a lot according to the 
chosen range of km'- in bright quasars (10-1 Jy), vm is mostly 
around lO^-lO-'GHz if the emitting regions are more com- 
pact and closer to the AGN core (case C2Co), and around 10- 
100 GHz in the other case (case C2Ex); fainter quasars have 
Vm ^ 400 GHz in the case C2Co and vm ^ 40 GHz in the case 
C2Ex. In general, the two cases give quite different results for 
the value of vm- For C2Co the majority of Jy and sub-Jy FSRQs 
can keep flat spectra up to v » 100 GHz; on the contrary, for 
C2Ex, apart from the few brightest quasars in the sky, the transi- 
tion to the optically thin synchrotron spectra can occur already at 
around 10-20 GHz. As a comparison, we also plot in Fig.Qthe 
break frequency value, vm - 60 GHz, of 3C 273 e stimated by 
multifrequency quasi-simultaneous observations in IClegg et aU 
(1983). This value is well within the range expected from the 
model C2Ex. Additional indications of a break in the spectra of 
bright blazars around 100 GHz have been re ported by analyses 
of Planck data (IPlanck CoUaborationil20TTillM) . 

Another physical quantity we can derive from the model is 
the magnetic field strength in the AGN jet, in equipartition con- 
ditions. We plot in Fig. |8] the values of //gq at the jet radius of 
1 pc for FSRQ and BL Lac sources with 5-GHz flux densities of 
55 = 1, and 0.1 Jy. Magnetic fields in BLLacs are found to be 
on average lower than in FSRQs, as expected in less powerful 
objects. The distributions for H^a, are quite sharp, with peak val- 
ues around 200-lOOOmG for FSRQs, and 100-500mG for BL 
Lacs. Magnetic field strengths of the order of few hundreds of 
mG are the values expected using the Konigl model in equipar- 
tition reg ime for sources o f J y or sub-Jy flux density (e.g., see 
results in lJiang et al.ll 19981 and lO^Sumvan & Gabuzdall2009l) . 
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5. High-frequency peaked spectrum sources 

The model described in Sect.|4]does not apply to sources with 
inverted spectra (a ^ 0.3). Very compact sources have con- 
vex spectra, peaking at GHz frequencies (GHz peaked spectrum, 
GPS, sources) or at tens of GHz (high frequency peakers, HFP). 
It is now widely agreed that GPS and HFP sources correspond 
to the early stages of evolution of powerful radio sources when 
the radio emitting region grows and expands in the interstellar 
medium of the host galaxy befo re becoming an extende d radio 
sourc e (see iDe Zotti et al.ll2005l and references therein; lO'Deal 
ll998l for a review). At low luminosity, inverted spectra may also 
correspond to late stages of AGN evolution, characte rized by low 
radiation/accretion efficiency (D e Zotti et a l.. 2010). 

About 5% of sources in the total NVSS/GB6 sample and 
~ 15% of flat-spectrum sources have inverted spectra between 1 
and 5 GHz; on the other hand, at high flux densities, these frac- 
tions increase up to about the 10% and 20% , respectively (see 
Table[TJ. The fraction of peaked sources increases significantly 
in higher freq u ency s urveys. For example, in the AT20G survey, 
iMurphv et all (1201 Ol) found that 21% of sources peak between 
5 and 20 GHz, and 14% have a spectrum rising up to 20 GHz 
(see also Hancock, 2009). If we restrict ourselves to the almost- 
complete subsample AT20G-dl5S100, we find about 22% of 
sources with Cj > 0.3, and 18% with > 0.3, of which about 
10% are inverted up to 20 GHz. 

Most of these sources cannot be interpreted as truly young 
GPS/HFPs, however Simultaneous multifrequency observations 
of high-frequency peaked or inverted-spectrum sources have 
given evidence of high flux-density variability or extended 
emission in a large number of them, which is not compatible 
with youth scena rio (Dallacasa et al., 2000; Stanahellini, 2003; 
Tomiainen etaO. [2005; Tinti et al., 2005; Bolton et al., 2006; 



Orienti et all 120061: iPlanck Collaborationii201 liL but indicative 
that they are most likely blazars caught in an active state, i.e. 
when a flaring, strongly self-absorbed synchrotron component 
dominates the emission spectrum. Moreover, when observations 
are not simultaneous, variability in flat-spectrum sources may 
lead to misinterpret them as inverted-spectrum sources. 

Our approach to this class of sources is purely statistical, 
without dealing with the physical origin of the emission. We 
are interested in the spectral shape of inverted-spectrum sources 
(i.e., all sources with a > 0.3 between 1 and 5 GHz) and in sta- 
tistically characterizing their behaviour at v > 5 GHz. 

In order to make predictions on the high-frequency number 
counts, we need some information on the peak frequency (Vpeak) 
and on the spectral index above the peak frequency. In the lit- 
erature we found three samples of sources with simultaneous 
measurements of a wide frequency range, allowing us to recover 
this information: 



- iKovalev et alj j 19991 hereafter K99) present a nearly simul- 
taneous observation of extragalactic radio sources at six fre- 
quencies (0.96, 2.30, 3.90, 7.70, 11.2, and 21.65 GHz) with 
the RATAN-600 radio tel escope. The sample c onsists of 546 
sources selected from the lPreston et al.l(ll985l) VLBI survey 
with a correlated flux density exceeding 0.1 Jy at 2.3 GHz. 
It includes primarily compact flat-spectrum objects because 
the Preston et al. sample is complete only for sources with 



-0.5. 



- iDaUacasa et"an (l2000l hereafter DOO) provide a list of 55 
high frequency peakers HFPs. This sample was obtained by 
selecting sources with flux density Sschz ^ 300 mJy and 

> 0.5 from the Green Bank survey 87GB and NVSS. 
Then HFP candidates were observed with the VLA at 1.4, 
1.7, 4.4, 5, 8, 8.5, 15, and 22.5 GHz, to acquire simultane- 
ous measurements of radio spectra and to remove variable 
sources that do not satisfy the criterion a > 0.5. 

- iBolton et alj ( |2004 hereafter B04) have carried out a follow- 
up of 176 sources from the 15-GHz 9C survey. Sources 
were selected from two complete samples with flux limits 
of 25 and 60 mJy. Simultaneous observations were made at 
frequencies of 1.4, 4.8, 22, and 43 GHz with the VLA and 
at 15 GHz with the Ryle Telescope. In addition, 51 sources 
were also observed at 3 1 GHz with the OVRO telescope. 



Spectra of sources from these three samples were first fit- 
ted with a single power law. If the fit was poor, we used 
a broken power law or a second-order polynomial in log- 



arithmic 



Then we selected all sources with a. 



space. „^ ^w,^^ ..lu. 

log(55/5i)/log(4.8/1.4) > 0.3 (whereas and5i are the values 
of the fit at 4.8 and 1.4 GHz). In this way, we created a sample 
of 142 inverted-spectrum sources: 78 of them come from K99, 
50 from DOO and 30 from B04 (16 are in common between K99 
and DOO). Apart from sources of B04, we were able to estimate 
the peak frequency only if it occurs at frequencies v < 20 GHz. 
We found that most of spectra peak at v ^ 10 GHz, and that 
117 sources (~ 80%) have Vpeak < 20 GHz. In the B04 sam- 
ple, where measurements extend up to 43 GHz, there are two 
objects with 20 < Vpeak < 43 GHz and another three that are 
still inverted at 43 GHz. In Fig.|9] we report the distribution for 
Vpeak up to 20 GHz fr om our sample and from the GPS/HFP can- 
didates collected by iTorniainen et al.l (l2008l) : the two distribu- 
tions agree well. In the plot we considered only the 41 sources 
in the Tomiainen et al. sample classified as "gps" or "convex" 
and with pea k frequency VpRA > 4 GHz. A similar distribution is 
also given bv lVollmer e'tanl2008h for a sample of 91 GPS/HFP 



This peak frequency, characteristic of GPS/HFP sources, must be 
clearly distinguished from the peak frequency, v,„ in the spectra of 
blazar sources, discussed in Sect. 4. 



candidates (see Fig. 7 in that paper). 

Results on the distribution of the peak frequency are sum- 
marized in Table[3] also for frequencies higher than 20 GHz. We 
recall that we select only sources with an inverted spectrum be- 
tween 1 and 5 GHz, and so the distribution is not meaningful 
for GPS sources with Vpeak < 5 GHz. Moreover, in Table[3] we 
assume that there are no inverted-spectrum sources with peak 
at frequencies higher than 100 GHz. In any case, current data 
indicate that the percentage of sources with Vpeak > 100 GHz 
has to be very small and is negligible for our purposes. This 
hypothesis is supported by the spectral behaviour of sources 
from simultaneous observations at 20 and 95 GHz by ATCA 
jSadler et all l200a) and by WMAP (iGonzalez-Nuevo et all 
20Ci). F inaUy, in simultaneous observations o f GPS/HFP candi- 
dates by ITorniainen et"an (I2005L l2007l 1200 8*) the maximum ob- 
served peak frequency is ~ 46 GHz, even if their measurements 
extend up to 250 GHz. 

For 98 objects of our sample of inverted-spectrum sources 
we can provide a reliable estimate of the spectral index after 
the peak (ff/,,). The spectral index distribution, plotted in Fig.|9] 
peaks at o-/,,- ^ -0.7, but it also has a lar ge tail of steep value s 
up to -2. Similar r esults were obtained by [Snellen et al. l (Il998l) . 
On the other hand, Ide Vries et akl (11997') found an average a/,,- ~ 
-0.7 but no sources with or;,, < -1 on a sample of 72 GPS. 
Since very few data are currently available on this subject, we 
use the spectral index distribution of Fig.|9]to predict the spectral 
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Fig. 9. Left panel: distiibution of the observed p eak frequencies for the sample of inverted-spectrum sources defined in the text 
(solid line), and for the GPS/HFP candidates in iTorniainen et alj (l2008l) . Right panel: distribution of the spectral indexes above 
the peak frequency for the inverted-spectrum sources. The solid line is given by f{a) - -A{a - ao) exp[-0.5(a - aof' jcr^^], with 
ao - -0.10 and ctq, - 0.53. 



Table 3. Distribution of the peak frequencies for sources with 
a\ > 0.3. 



Peak Frequency (GHz) 


Source Fraction 


Vpeak < 6 


0.27 


6 < Vpeak < 12 


0.40 


12 < Vpeak < 20 


0.15 


20 < Vpeak < 40 


0.15 


40 < Vpeak < 100 


0.03 



behaviour of inverted-spectrum sources (a > 0.3) at frequencies 

above Vpeak- 

6. Extrapolation of 5-GHz number counts to higher 
frequencies 

Based on the 5 GHz number counts and on information of spec- 
tral properties of ERS described in previous sections, we are now 
able to make predictions of number counts at cm/mm wave- 
lengths. We deal with the complexity of source spectra as fol- 
lows. 

- We consider three different populations of radio sources, ac- 
cording to their spectral index in the frequency range 1- 
5 GHz: steep-spectrum if ffj < -0.5; flat-spectrum if -0.5 < 
OTj < 0.3; inverted-spectrum if aj > 0.3. Two flat-spectrum 
populations are taken into account: FSRQs and BLLacs. 

- A simulated source catalogue is produced at 5 GHz. For flat 
plus inverted-spectrum sources we use the n(S) obtained 
from the fit of observational data given by Eq. ([T]l. The n{S) 
of steep-spectrum sources is then extracted a s the difference 
between the total number counts from the iToffolatti et al.l 
([1998) model and the fit of Eq. O; see also Fig.|2] The 
fraction of BL Lacs among flat-spectrum sources at 5 GHz, 
which depends on the flux density, is taken from the evolu- 
tionary model of De Zotti et al. (2005). 

The spectrum for each source is approximated as a power law 
at V < 5 GHz. Spectral indices a^. are supposed to follow the 



flux-density dependent spectral index distributions obtained 
from the NVSS/GB6 sample (see SectionlTTland Tabled- 
The 5-GHz flux density of each source is then extrapolated 
to V > 5 GHz using the following spectral models. 

For steep-spectrum sources we take 



S(v) = S5cm(v/5r 



(9) 



where a/,, - ffj -Aa. The spectral steepening Ao- is randomly 
extracted from a Gaussian distribution with average 0.3 and 
dispersion 0.2, in agreement with the steepening found in the 
AT20G-dl5S50 sample and shown in Fig.|4] A small per- 
centage of flattening or upturning sources is also included. 

For inverted-spectrum sources we use 



5(v) = 5o(y/vo)'(l-e-<''/^°'"), 



(10) 



where vq is directly related to the peak frequency in the spec- 
trum and Sq - (1 -e"')5(vo). The coefficients k and / are the 
spectral indices for the rising and declining parts of the spec- 
trum. We set A: = ffj, whereas / is extracted from the best-fit 
distribution shown in Fig.|9] The peak frequency is extracted 
from the distribution in Table|3] 

For flat-spectrum sources we use different spectral models. 
In the simplest one (CO) radio spectra at v > 5 GHz are mod- 



Aa is extracted from Gaussian distributions with 



elled as power laws with spectral index a/i - a; - Aa, where 



< Aa >- 0.1 and cr - OA 

< Aa >= 0.2 and cr = 0.2 



if -0.5 < 0-= < -0.1 



if 



•0.1<a;<0.3. (11) 



These values for < Aa > (and the associated disper- 
sions) agree with observations from the AT20G survey (see 
Section[3l1i. 

As another step in modelling the spectra of flat-spectrum 
ERS, we introduce a "break frequency", as discussed in 
Sec.m Spectra are now described by 



Sscnziv/Sr" 
S(vm) (v/vmT" 



5 < V <vm 

V>Vm, 



(12) 
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Fig. 10. Predicted differential number counts normalized to S~^'^ at 8.44 GHz (Left panel) and at 15 GHz (Right panel) for the 
CI model (thick continuous line), i.e. the model with no difference between t he value adopted f or FSRQs and BLLac objects 
(see text). At 8.44GHz data points are a collection from different samples (see lDe Zotti et al.ll2010l) . whereas at 15 GHz they have 
been computed from the IOC survey. In the left panel predictions for number counts of flat- (short-dashed line) and steep-spectrum 
(dotted line) sources are also shown, compared with results from measurements at 8.4 GHz of flat-spectrum sources by the CRATES 
program (empty circular points). Predictions at 15 GHz from the lDe Zotti et al.l ( I2005h model are also shown (blue dash-dotted hne). 



where vm is the break frequency, which is interpreted as the 
frequency at which the transition from the optically thick to 
the optically thin regime of synchrotron radiation occurs. As 
above, the spectral index in the optically thick regime at v > 
5 GHz is a// - a^^-Aa, with Aa extracted from the Gaussian 
distributions of Eq. (fTTl i. whereas a,, is the spectral index 
in the optically thin regime, determined by synchrotron or 
IC energy loss es, and is taken to be a,, - -0.8 + 0.2, in 
agreement with iKellermannI (1 1 9661) . 

The break frequency is related to the value of the parameter 
rM, discussed in Sec. l4.1l We considered three cases: 

- case CI: 0.01 < < 10 pc, with no difference between 
FSRQs and BL Lac objects; 

- case C2Co: 0.01 < < 0.3 pc for BL Lacs and 

0. 03 < < 1 pc for FSRQs, i.e. the radio emission 
in FSRQs arises from a "more compact" region than in 
the following case, C2Ex; 

- case C2Ex: 0.01 < tm < 0.3 pc, the same value as before 
for BL Lac objects, and 0.3 < tm < 10 pc for FSRQs, 

1. e. the radio emission in FSRQs comes from a relatively 
"more extended" region than in the case C2Co. 

- In the end, simulated catalogues can be extracted at different 
frequencies and for different flux limits. Number counts and 
spectral properties of the three populations of ERS can then 
be compared with observational data. 

7. Predictions on number counts and spectral 
properties of ERS at v > 5 GHz 

7.1. Number counts: predictions vs observations 

A summary of the relevant data is presented in AppendixlAl 
At V ^ 30 GHz, we find a general agreement between the ob- 



served n(S ) and predictions of all models described in the pre- 
vious section (CO, CI, C2Co or C2Ex). In Fig.[TO]we plot our 
results at 8.4 and 15 GHz, but only for our intermediate model, 
i.e. CI, because of the very small differences among models. In 
Fig.[TT] thanks to the almost-simultaneous measurements at 5 
and 20 GHz for most of the AT20G sources, we test in more de- 
tail the goodness of our predictions for the total number counts 
of sources in the sub-sample AT20G-dl5S100 (left panel), as 
well as for sources within different ranges of af^ (right panel). 
Again, all our models are found to agree very well with the 
data. This is an important confirmation that the adopted num- 
ber counts at 5 GHz (especially for flat-spectrum sources) and 
spectral index distributions from the NVSS/GB6 sample are es- 
sentially correct. 

The agreement with the data is confirmed at 30-33 GHz over 
a very broad flux density range (Fig. [12}, and it extends up to 
100 GHz, as shown in Fig.[T3] 

To distinguish among our models we need data at still higher 
frequencies (> 100 GHz; Fig.[T4li. The model without a spectral 
break (CO) overestimates the counts over essentially the whole 
flux-density range. The other models show the largest differ- 
ences among each other at the brightest fluxes (at S = 1 Jy, the 
number counts from the model C2Co are almost a factor 2 higher 
compared to the model C2Ex), whereas they tend to converge at 
S < lOmJy. A very good match with the data is obtained with 
the model C2Ex, whereas the other models clearly overestimate 
the Planck counts. 

For a direct compa rison, we also plot in Fig.fT4lthe number 
counts predicted by the lDe Zotti et al.l (|2005|) cosmological evo- 
lution model for ERS. At bright flux densities (S > 0.1 Jy), their 
results are compatible with our model CO, as expected because 
both models use simple power-law approximations for blazar 
spectra. However, at lower fluxes, the differences between these 
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Fig. 11. Normalized differential number counts at 20 GHz from AT20G-dl5S 100 sub-sample of the AT20G survey (filled squares; 
see text) compared to predictions of the models described in the text: C O(dotted lines), CI (th ick continuous lines), C2Ex (lower 
red long-dashed lines) and C2Co (upper red long-dashed lines) and the iDe Zotti et alj (l2005l) model (blue dash-dotted line). Left 
panel: the differential number counts for all source populations in the sample. The empty circles are from the 23-GHz WMAP 
sources catalogue. Right panel: differential number counts for AT20G-dl5S100 sources with -0.5 < of^" < 0.3 (filled squares), 
a^^ < -0.5 (empty squares) and o-^" > 0.3 (star points) and the corresponding model predictions. 



two models increase, and the De Zotti et al. model provides pre- 
dictions similar to those obtained from our models including a 
break frequency. This is no surprises because it can be attributed 
to the different 5-GHz n{S ) for flat-spectrum sources that our 
models use (as shown in Fig.|2|i, if compared to the De Zotti et 
al. model 

Finally, a clear distinction among our models should be pro- 
vided by the radio source counts at the highest Planck HFI fre- 
quencies. Our predictions on integral number counts, A^(> S), 
are presented in Fig.[T5]and in Table|4] As an example, the num- 
ber of objects brighter than 1 Jy should be reduced by ~ 40% 
in the range 353-857 GHz for the model C2Ex, but by only 
~ 20% or less in the other cases. Moreover, we predict that 
the Herschel-AILAS survey (Ba les et all l2010) will detect from 
25 (model C2Ex) to 50 (model C2Co) blazars brighter than 
50mJy at 500 um ove r its area of 550deg^. For comparison, the 
( De Zotti et al. , 2005b model y ields ~ 80 blazars in the same area 
(iGonzalez-N uevo et al.Ll2010l) . 

7.2. Predicted spectral properties of radio sources 

As a further test, we considered the average and the median 
spectral indices of ERS in catalogues selected at high frequen- 
cies. In the ACT sample, all but two of the 42 detections with 
flux density ^ug > 50 mJy have counterparts in AT20G at 5 
and 20 GHz (Marr iage et ail 1201 ih . An average spectral steep- 
ening is found above 20 GHz: the median spectral index and 
the dispersion of the distribution are a^^ - -0.07 + 0.37, 




0.001 
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Fig. 12. Predicted differential number counts at 33 GHz com- 
pared to observational data. The lines have the same mean- 
ing as in Fig.[TT] The area labelled 'GBT' shows estimates by 
■Mason et al. (2009), based on the 31-GHz survey conducted us- 
ing the GBT and OVRO telescopes. 



' The excess of steep-spectrum ERS - with respect to the observed 
ones - adopted by the De Zotti et al. model at 5 GHz can clearly mimic 
the effect of the break frequency in the spectra of flat-spectrum ERS at 
mm wavelengths at least partially. 



,,148 



-0.20 : 



:0.21 and 0-148 



-0.39 + 0.24. Somewhat flatter 



spectral indices are observed for SPT sources (57 objects with 
a 5-cr detection at 148 and 220 GHz) between 5 and 148 GHz, 
< Qf^'18 >= -0.13 + 0.21, with a substantial steepening between 
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100 




B 1 



0.1 




100 GHz 



WMAP 

Planck 




S(Jy) 



Table 6. Predicted median spectral indices of ERS between 
30 GHz and the Planck frequency channel indicated in the first 
row. The corresponding median values calculated for the spec- 
tral indices of ERS in the P lanck ERCSC are taken from the 
iPlanck Collaboration! (1201 lib . 



v[GHz] 


44 


70 


100 


143 


217 


Planck 


-0.06 


-0.18 


-0.28 


-0.39 


-0.37 


model 












CO 


-0.11 


-0.13 


-0.14 


-0.14 


-0.14 


CI 


-0.18 


-0.23 


-0.27 


-0.31 


-0.34 


C2Co 


-0.13 


-0.17 


-0.19 


-0.21 


-0.24 


C2Ex 


-0.22 


-0.28 


-0.34 


-0.39 


-0.44 



V ^ 100 GHz) give quite "flat" mean 5-148 GHz spectral indices 
(<_a\'^^ >- -0.2), very close to the value found by Vieir a et aH 
( 2010). This is because of the 148 GHz selection that emphasizes 
sources with flatter spectra. 

On the other hand, the median spectral indices yielded by 
models CO and C2Co are not compatible with observations at 
frequencies ^ 100 GHz. 



Fig. 13. Differential number counts predictions as in Fig.lTTIbut 
at 44, 70 and 100 GHz Bigger black points a re from the WMAP 
NEWPS catalogue bv lMassardi et all (l2009l) at 41 and 61 GHz; 
smaller blue points are from the Planck ERCSC. For clarity, 
number counts at 70 GHz are multiplied by 0.1, at 100 GHz by 
0.01. 

Table 4. Number of ERS brighter than 1 and 0.5 Jy in the fuU 
sky at the highest Planck frequencies predicted by our models. 







A'(> 1 Jy) 




A'(> 0.5 Jy) 


v[GHz] 


353 


545 


857 


353 


545 857 


model 












CO 


205 


207 


218 


507 


512 522 


CI 


130 


116 


107 


311 


278 246 


C2Co 


174 


156 


141 


406 


361 313 


C2Ex 


77 


60 


47 


184 


143 112 



148 and 220 GHz, < a^zo >^ _o.5 dVieira et al.Ll2010l) . To com- 
pare these data with model predictions, we produced simulated 
catalogues of sources selected at 148 GHz with the same flux 
limits as the SPT and ACT catalogues and with S 20 > 50 mJy, 
and we computed the mean and median spectral indices among 
pairs of frequencies in the range 5-220 GHz; results are reported 
in Table|5l 

Similarly, we used the models to simulate samples selected 
at 30 GHz with a completeness limit of 1 Jy to compare the pre- 
dicted median values of spectral indices between 30 GHz and 
higher Planck frequencies and the associ ated dispersions with 
the obs erved values, given in Table 2 of the lPlanck Collaboration! 
(1201 lib . The results are shown in Table|6] 

In general, we see that models including a break frequency 
in blazars spectra are required to explain the median spectral 
indices ^ -0.3 observed when the surveys are performed at fre- 
quencies > 100 GHz. In particular, the models CI and C2Ex 
yield median spectral indices very similar to the Planck ERCSC 
ones, i.e. -0.4 ^ $ -0.3. They also yield median steep- 
enings af - a'^s 9.14 (CI) and 0.26 (C2Ex) to be com- 
pared with 0.32 from ACT data, and a^'^^ - of 0.22 and 
0.28, respectively, to be compared with 0.37 from SPT data 
(see Table|5]l. Remarkably, models including a spectral break (at 



8. Conclusions 

The main goal of this paper was to present physically grounded 
models to extrapolate the number counts of ERS, observationally 
determined over very large flux density intervals at cm wave- 
lengths down to mm wavelengths, where the majority of the ex- 
periments aimed at accurately measuring CMB anisotropies are 
carried out. Accurate extrapolations are necessary to minimize 
and/or to control the contamination of CMB maps by unresolved 
extragalactic sources. Moreover, this paper makes a first attempt 
at constraining the most relevant physical parameters that char- 
acterize the emission of blazar sources by using the number 
counts and the spectral properties of ERS estimated from high- 
frequency radio surveys. 

We focussed on spectra of blazars that dominate the 
mm-wave number counts of ERS at bright flux densities. 
In the most recent data s ets (e.g., Massardi et al., 201()3 
iPlanck Collaboration! |20I lilkh . a relevant steepening in blazar 
spectra with emerging spectral indices in the interval between 
-0.5 and -1.2, is commonly observed. We interpreted this 
spectral behavior as caused, at least partially, by the transi- 
tion from the optically-thick to the optically-thin regime in 
the observed sy nchrotron emission of AGN jets (see, e.g., 
lMarscheil[T996l) . Indeed, a "break" in the synchrotron spectrum 
of blazars above which the spectrum becomes "'steep" is pre- 
dicted by models of synchrotron emission from inhomogeneous 
unresolved relativistic jets ( Bl andford & K onigl. 1979; KoniglL 
1 198 lb iMarscher & Geal [ 19851) . We estimated the value of the 
frequency vm at which the break occurs on the basis of the ERS 
flux densities measured at 5 GHz and of the most typical values 
for the relevant physical parameters of AGNs. In the framework 
of these models, the most relevant and critical physical param- 
eter is the dimension of the region (approximated as homoge- 
neous and spherical) that is mainly responsible of the emission 
at the break frequency. For a conical jet model, this parameter 
can be related to the distance of the emitting region from the 
AGN core (r^; see Appendix A). 

We have investigated four possible cases. The simpler and 
basic one, CO, assumes simple power-law spectra above 5 GHz, 
with a dispersion of spectral indices. The other three models 
feature at high frequencies a spectral break and make differ- 
ent choices for the parameter r^. In model CI the range of 
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Table 5. Means, dispersions, and medians of spectral indices calculated between different pairs of frequencies as predicted by our 
models. 



y[GHz] 


[5,20] 


[20, 148] 


[5, 148] 


[5, 148] 


[148,220] 






ACT (median ±cr) 




SPT (mean ±cr) 




-0.07 ± 0.37 


-0.39 ± 0.24 - 


-V.zu ± U.Zl 


-0.13 ±0.21 


-0.50 


model 


ACT simulated sampl 


e 


SPT simulated sample 


< a > ±aa CO 


-0.08 ± 0.30 


-0.16 ±0.27 - 


-0.13 ±0.26 


-0.08 ± 0.29 


-0.14 ±0.32 


CI 


-0.08 ± 0.33 


-0.25 ± 0.30 - 


-0.18 ±0.28 


-0.16 ±0.31 


-0.37 ± 0.38 


C2Co 


-0.07 ±0.31 


-0.20 ± 0.28 - 


-0.15 ±0.27 


-0.14 ±0.30 


-0.36 ± 0.39 


C2Ex 


-0.07 + 0.36 


-0.35 ± 0.31 


0.23 ± 0.30 


-0.22 ± 0.33 


-0.50 ± 0.37 


a (median) CO 


-0.10 


-0.17 


-0.13 


-0.09 


-0.15 


CI 


-0.10 


-0.24 


-0.18 


-0.17 


-0.37 


C2Co 


-0.09 


-0.21 


-0.15 


-0.15 


-0.37 


C2Ex 


-0.09 


-0.35 


-0.24 


-0.24 


-0.55 




Fig. 14. Comparison between predicted and observed differential number counts at 148 GHz (left panel) and at 220 GHz (ri^ht 
panel) . Filled circles: ACT data; open black circles: SPT data; open blue circles: Planck ERCSC counts (iPlanck CollaborationL 
l2011ih at 143 GHz (left panel) and 217 GHz (right panel). The lines have the same meaning as in Fig.fTTl 



rjvf (and therefore of vm) is very broad and should include all 
most likely valu es expected from theoretical models (see, e.g. 
iLobanovl |2010l for a recent review). Models C2Co and C2Ex 
assume different distributions of rM for BLLacs and FSRQs, 
with the former objects that generate, in general, the synchrotron 
emission from more compact regions, implying higher values of 
Vm (above 100 GHz for bright objects). These two models differ 
only in the r^j distributions for FSRQs: in the C2Co model the 
emitting regions are more compact, implying values of vm par- 
tially overlapping with those for BLLacs, whereas in the C2Ex 
model they are more extended and the values of vm are con- 
sequently lower and clearly distinguishable from those adopted 
for BLLacs (approximately in the range of 10-300 GHz for 
S > 1 Jy). 

Obviously, high frequency (v > 100 GHz) data are the most 
powerful for distinguishing among those models. These data 
clearly require spectral breaks, thus ruling out model CO, in spite 
of the average steepening that was introduced in the spectral in- 
dices of ERS at V > 5 GHz, and they favour the model C2Ex. 
According to this, most of the FSRQs (which are the dominant 
population at low frequencies and at Jy flux densities), differ- 



ently from BLLacs, should bend their flat spectrum before or 
around 100 GHz. The model also predicts a substantial increase 
of the BL Lac fraction at high frequencies and bright flux densi- 
ties. This is indeed observed: a clear dichotomy between FSRQs 
and BL Lac objects has been found in the Planck ERCSC, where 
almost all radio sources show very flat spectral indices at LFI fre- 
quencies, i.e. aifj ^ -0.2, whereas at HFI frequencies, BL Lacs 
keep flat spectra, i.e. uhfi ~ -0.5, and a high fraction of FSRQs 
show steeper spectra, i.e. anpi < -0.5. Moreover, the fraction of 
BL Lacs above 1 Jy increases from 10% in the 5-GHz lKuhr et al.l 
( 198lb sample t o 20% in the 37 GHz selected sample discussed 
by the iPlanck Collaboratio ^ (|20TT1 . More constraints are ex- 
pected f rom sub-mm counts of blazars produced by the Herschel 
surveys (iGonzalez-Nuevo et al I I20T0I) . 

According to our model, these results imply that the pa- 
rameter rM should be of parsec-scales, at least for FSRQs, 
in agreement with theoretical predictions (iMarscher & Gear! 
1 985) , whereas values of <K 1 pc should be only typical of 
BLLac objects or of rare quasar sources. 

On the other hand, the model C2Ex slightly underestimates 
the number counts at 20-30 GHz and predicts too steep median 
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S(Jy) 



Fig. 15. Predictions on integral number counts at 353, 545 and 
857 GHz from the models CO, CI, C2Co, C2Ex discussed in 
Section 7. The plotted lines have the same meaning as in Fig.fTTI 



spectra at the Planck LFI frequencies. This indicates that spec- 
tral breaks at v < 20 GHz are rarer, i.e. the emitting regions 
of FSRQs, particularly of the fainter ones, are somewhat more 
compact than implied by this model. The possibility of more 
compact AGN cores associated with less powerful sources could 
agree better with physical models of the AGN jet emission (e.g., 
[Ghisellini et al., 1998), and could provide an improvement of the 
model predictions at tens of GHz. 

The physical model used in this paper adopts only a simpli- 
fied description of the synchrotron emission in AGN jets and, of 
course, does not pretend to take into account the complexity of 
all physical mechanisms involved. Nevertheless, it is capable of 
providing a clear interpretation of several features observed in 
blazar spectra. The very interesting quantitative results obtained 
on number counts and on spectral index distributions of ERS in 
the whole frequency interval from 5 to 220 GHz clearly support 
this conclusion. A step forward will be the re-analysis of the 
luminosity functions of blazar sources at radio frequencies, to 
follow the cosmological evolution of this class of sources for im- 
proving the fit, and the corresponding interpretation of the source 
number counts at v > 100 GHz and at faint flux levels. However, 
this further step will be the subject of a following paper. 
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Appendix A: Large-area surveys of radio sources 

A.1. Low-frequeny data: the catalogues of steep- and 
flat-spectrum sources 

Large-area deep surveys are specially present at GHz frequen- 
cies and have allowed us to make statistical studies of the spec- 
tral behaviour of ERS for the basic classification into steep- and 
flat-spectrum sources. 

- By us ing flie NRAO VLA Sky Survey (NVSS: ICondon et al.L 
19981) at 1.4 GHz and the Green Bank survey (GB6; 
Gregory et a il 1996) at 4.85 GHz, we built a large and 
complete catalogue of steep- and flat-spectrum sources at 

5 GHz. These two surveys overlap at declinations of 0° < 

6 < -75° (i.e., an area Q ^ 6.07 sr) and have allowed us 
to calculate the spectral indices of ERS in almost half of the 
sky. The resolution of the two surveys is quite different, be- 
ing 45 arcsec in NVSS and 3.5 arcmin in GB6, whereas the 
flux limits are 2.5 and 18 mJy respectively. 

At first, we decided to consider only GB6 sources with flux 
density 5 4.85 > lOOmJy and Galactic latitude \b\ > 10°. 
The 100-mJy flux limit for GB6 sources is well above the 
flux Umit of the two surveys and guarantees that we are not 
losing sources with a rising spectrum, a > Q. Moreover, it 
minimizes the effects of flux density errors. The Galactic cut 
at latitude < 10° is to avoid Galactic sources in the sam- 
ple. In addition, we excluded GB6 sources flagged as "W" 
(not reliable source) and "C" (confused source). We cross- 
matched the GB6 source positions with the NVSS catalogue 
by taking all positive matches within a position offset of 
89 arcsec. We used a shghtly bigger maximum position offset 
with respect to lHealev et al.l ([2007 ) for reducing the number 
of GB6 sources without NVSS counterparts. In this way, we 
found 8127 sources with a counterpart in NVSS and only 
23 without. Finally, whenever more than one of the NVSS 
sources fell within the GB6 beam, as a consequence of the 
better NVSS angular resolution, which may lead to individu- 
ally resolved multiple components, we summed their fluxes, 
correcting for the effect of the GB6 beam. We ended up 
with 2975 flat-spectrum sources (corresponding to ~ 37% of 
the total number of sources in the sample) and 5152 steep- 
spectrum sources (~ 63%). Hereafter, we indicate this source 
catalogue with spectral information as our NVSS/GB6 sam- 
ple. 

- The CRATES programme (iHealev et al.L l2007l) has also car- 
ried out an almost all-sky sample of flat-spectrum sources 
brighter than 65 mJy at 5 GHz, using the existing surveys 
at GHz frequencies. Then they assembled the 8.4-GHz flux 
densities of flat -spect r um so urces from observations done by 
CLASS (■Mverserani200l and from new observations by 
VLA and ATCA. To mantain uniformity in our analysis, we 
took into account only CRATES sources in the GB6 area, i.e. 
where spectra are computed from NVSS and GB6 measure- 
ments. In this area, the authors found ~ 5000 flat-spectrum 
sources with 5-GHz flux density 5 > 65 mJy (about a 33% 
of total sources). 

- A sample of flat-spectrum s ources is also pr ovided by 
the Parkes quarter-J y sample (I Jackson et al.l 12002 : see also 
lDeZottiet"al]|2005h . It consists of 878 objects selected at 
2.7 GHz from several complete surveys of the Parkes radio 
source catalogue, and having spectral index between 2.7 and 
5 GHz Q'2 7 > -0.4. The flux limit of these surveys varies be- 
tween 0. 1 and 0.6 Jy, although it is 0.25 Jy for most of them. 
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A.2. Surveys at frequencies higher than 10 GHz 

Recent experiments have sui-veyed large areas of the sky at fre- 
quencies higher than 10 GHz. They are important to test the va- 
lidity of our predictions on number counts of ERS, but also to 
provide more direct information on the spectral shape of sources 
at cm/mm wavelengths. 

- The Ryle-Telescope 9C surveys (iTavlor et al.L l200lt 
IWaldram et al.L l2003h have provided a catalogue of sources 
at 15 GHz with a completeness limit of 25mJy. These sur- 
veys cover the fields observed by the Very Small Array 
(VSA), corresponding to an area of ~ 520 deg^. Moreover, 
IWaldram et al.l (|2009|) have reported on a series of deeper 
regions, amounting to an area of 115 deg^ complete to ap- 
proximately lOmJy, and of 29deg^ complete to approxi- 
mately 5.5 mJy. Finally , the Tenth Cambridge (IOC) Survey 
{ami Co nsortium! |20 1 0) has covered an area of » 27 deg^ at 
15.7 GHz down to a completeness limit of 1 mJy (within it, 
some deeper areas, covering x 12deg^, are complete down 
to 0.5 mJy). 

- A 20-GHz survey of the whole southern sky has been car- 
ried out by the Australian Telescope Compact Array (ATCA) 
from 2004 to 2008. The fiiU sour ce catal ogue (AT20G) 
is presented in iMurphv et al.l (1201 Oh and in iMassardi et al.l 
f2010b), and it includes 5890 sources above a flux-density 
limit of 40 mJy. The completeness of the AT20G catalogue 
is 93 per cent above 100 mJy and 78 per cent above 50mJy 
in regions south of declination -15°. Most of sources with 
declination 6 < -15° were also followed-up at 5 and 8 GHz 
by near-simultaneous observations. 

In our analysis we made an extensive use of this survey 
but we limited ourselves to the almost-complete sample of 
the AT20G catalogue at declinations 6 < -15°. It consists 
of 2195 sources with flux density S > lOOmJy and of 
1612 with 50 < 5 < 100 mJy. We called these two sub- 
samples AT20G-dl5S100 and AT20G-dl5S50, if the flux 
limit is lOOmJy and 50mJy respectively. For these sources, 
we were able to determine the spectral index at low frequen- 
cies by exploiting ATCA 5-GHz measurements and low- 
frequency surveys in the southern sky: the NVSS at 1.4 GHz 
for 6 > - 40°; the Sydney Un iversity Molonglo Sky Survey 
(SUMSS iMauch at aLL l2003h and the Molonglo Galactic 
Plane Survev (MGPS Murp hv etalL |2007) at 843 MHz for 
6 < -40°. Whereas 5-GHz ATCA measurements were not 
available, we searched for the source count erparts in the 
Southern Pa rkes-MIT-NRAO (PMN) survey dWright et al.L 
11994 Il996h . In AT20G-dl5S100 we obtained the spectral 
index for 2158 sources, more than 98% of the sub-sample: 
there were 395 sources with steep spectiTim and 1763 with 
flat spectrum (~ 82%), of which 471 with inverted spectiTim 
a > 0.3. In Table lA.ll we report the number of sources with 
5-GHz measurements and the number of steep- and flat- 
spectrum sources in the two sub-samples of the AT20G. The 
percentages of fl at- and steep-spectrum sources agree with 
the ones given bv lMassardi et al.l (l2010bh over a slightly dif- 
ferent and smaller area of AT20G (see their Figure 1), where 
they found ~ 82% and ~ 74% of flat-spectrum sources with 
flux density > 100 and 50mJy respectively. 

- At V ^ 30 GHz various CMB experiments have provided 
samples of extragalactic radio sources, giving estimates of 
the number counts at different ranges of flux densities: the 
Cosmic Background Imager (CBI; Mason et al., 2003) in the 
range 5-50 mJy; the Degr ee Angular Scale Interferometer 
(DASL lKovac et alil2002h for 5 > 100 mJy; the Very Smafl 



Table A.l. Spectral information in the two almost-complete 
sub-samples of AT20G. Columns give the number of sources; 
number of sources with estimated 1-5 GHz spectral index; num- 
ber of flat- and steep-spectrum sources. 



Sample 




N(a]) 


flat 


steep 


AT20G-dl5S100 


2195 


2158 


1763 


395 


AT20G-dl5S50 


3807 


3699 


2821 


878 



Array (VSA; IClearv et al.l 120051) in the range 20-114nJy ; 
the Sunyaev-Zel'dovich Array (SZA: lMuchovei et al.Ll2010h 
in the range 0.7-15 mJy. 

The Wilkinson Microwave Anisotropy Probe (WMAP) car- 
ried out all-sky surveys at 23, 33, 41, 61, and 94 GHz and 
provided a catalogue of ERS at a completeness levels of 
<: 1 Jy. Analyses from the WMAP team have yielded 390 
point sources in the five-year data (Wiight et al., 2009|), 
wher eas 62 new point sources are found in the seven-year 
data dGoldetaLLlToiO). WMAP five ^ear maps has been 
also analyzed by Massardi et al.l (l2009h that detected 516 
point sources, 457 of which were previously identified as ex- 
tragalactic sources. 

Data on ERS are also present fo r frequencies v ^ 1 00 GHz: 
the South Pole Telescope (SPT; IVieira et al.L 1201 Oh carried 
out a survey of ERS at 1.4 and 2.0 mm wavelengths with 
arcmin resolution and mJy depth over an area of87deg2;the 
Atacama Cosmology Telescope (ACT; Marriage et al., 201 1]) 
provided a catalog of 157 sources with flux density between 
15 and 1500 mJy det ected at 148 GHz in an area of 455 deg^. 
The Planck ERCSC dPlanck CoUaborationi |20TT3) reported 
data on compact sources detected in the nine Planck fre- 
quency channels between 30 and 857 GHz during the first 
1 .6 full-sky surveys. The analysis of the P lanck ERCSC data 
presented in lPlanck Collaboration! (1201 l i*) is limited to a pri- 
mary sample of 533 compact extragalactic sources at \b\ > 
5°, selected at 30 GHz. More than the 97% of these com- 
pact objects have been identified in external, published cat- 
alogues of ERS at GHz frequencies (see the Planck ERCSC 
Explanatory Supplement, for more details). Moreover, this 
30-GHz sample is found to be statistically complete down 
to a flux density of « 0.9-1.0 Jy and 290 ERS are found at 
above this flux density limit. The Planck Collaboration has 
been able to measure the number counts at the Planck LFI 
(30, 44 and 70 GHz) and at the HFI frequencies of 100, 143 
and 217 GHz, with an estimated com pleteness limits of 1.0, 
1.5, 1 .1, 0.9, 0.5, 0.4 Jy respectively (iPlanck CoUaborationi 
I2011ih . 



Appendix B: Estimate of the break frequency in 
blazars 

Spectra for the synchrotron emission from a spherical and ho- 
mogeneous source have a peak due to the self-absorption of 
their own radiation. The observed frequency at which the peak 
of the self-absorption occurs depends on the magnetic field and 
the de pth of the source, and it can be computed by dPacholczvfl 
U970h 



^syn, abs 



(B.l) 

where 6 is the observed angular dimension of the source and 
Vsyn,abs IS measured in GHz. The parameter p depends on the 
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emission model and gives the enhancement of the observed flux 
due to the beaming (Sobs - 6''^'^ S). It is equal to 3 for a moving 
isotro pic source and 2 for a cont inuous jet (see Ghisellini et al. 
ll993l and lUrrv & Padovanill 19951 for detailed discussion). In the 
analysis we used as reference value p = 3. We also assumed 
that emitting electrons have a power-law energy distribution 
N(y) = Ky-^'^-^"\ where a is the spectral index of the optically 
thin synchrotron emission. The term C„ in Eq. lB.ll depends on a 
by 

4c6(a) '^'^ 



where 



Ca - 2C] 



Tja)c5(a) 



(B.2) 



with ci = 3e/4mn c and t„, the optical depth of t he source at 
v,n. Th e functions csia) and c^ia) are provided in jPacholczvkl 

(ii97(y). 

As discussed in the text, the break frequency (vm) in a 
flat-spectrum source is approximately the synchrotron self- 
absorption frequency Vsyn,abs for the innermost part of the jet 
whose emission is observed at cm/mm wavelengths. If this re- 
gion is assumed to be homogeneous and spherical (with diameter 
d), Vm can be obtained from Eq. lB.ll For a conical jet geometry, 
the diameter is related to the distance from the AGN core by 
d - 2rM tan(0/2) ^ (prM, where (p is the semiangle of the conical 
jet. 

If the flux density for a flat-spectrum source is known at the 
observational frequency Vg, the observed flux density at vm (S m) 
can be extrapolated from using a power law spectrum (we are 
assuming that 5 m ^ S sy„,abs, i-e. the contributions from other jet 
regions are negligible at the frequency vm)- 

Sm^Cs with Cs = y7" S (vo) ■ (B.3) 

Moreover, we have seen from Eq.|5]Mlthat the magnetic field 
in equipartition condition is 



H,g ^ 1.06 X 10" 



(L/Watt) 



L(y/Kpc^)J 



2/7 



[mG] 



(B.4) 



The total luminosity of the source L can be calculated by the 
integral of the observed flux density: 



= J L(v)dv = :^6' 



(l+z) 



I+Q- 



S(v)dv [Watt], (B.5) 



where Dl is the luminosity distance in Mpc, Vmin and Vcut give the 
frequency range where the source emission is concentrated (we 
take v„,i„ - 10 MHz and v„„ = 10^ GHz, in agreement with stan- 
dard assumptions on synchtrotron emission in blazar sources). 
We used the relation between the luminosity emitted at a given 
frequency and the observed flux density, which takes into ac- 
count the K-correction and the relativistic beaming effects on 
the flux. The integral on the flux density is now expressed as a 
function of 5 m and vm' 



I S{v)dv = I S m\ — I dv 

Jv..„„ Jv„„„ ^VmI 

_ IVmm\ 
\ Vcut ) 



1 + a 



(B.6) 



for or -1 [if Of = -1 the integral is equal to S mVm ln(v„„/v,„,„). 
Below we assume this condition is verified; it is easy to extend 
the calculation for a - -1]. Using Eq. lB.3l the total luminosity 
becomes 



£)2ga-p 



[Watt], 



(B.7) 



Cl ^ 9.5 X W\nS{vo)Vo 



-ff// ' cut 



\ + a 



(Vmin \ 
Vein ) 



(B.8) 



By assuming for s implicity that the source is spherical w ith di- 
ameter ^ O.lrM (;Konig]lll98l]:lGhisellini & Tavecchioll2009h . 
we obtain the expression for the magnetic field: 



where 



-6/7 2(.,,-.)/7 



Ch ^ 0.47 X lO^^C^^^ . 



(B.9) 



(B.IO) 



The observed angular dimension of the source in Eg. IB. H is 



= C9(1+z>m/£)l [mas]. 



(B.ll) 



with Ce ^ 2.1 X 10^ 

FinaUy, including Eq. lB3tiR9tiBTT] in Eq.lO the break 
frequency becomes 



Vm = C(a, an) 



VP 



(B.12) 



where C{a, an) = (C„ ^Cg'^^^C]/^)'/^, and ^ = I + ^a 

16 „ o _ 32 ft _ 2g-51 r, _ 34fl_i,2„ 16 „ 

■^af 1,130 - 35, Pz - -^,Pd - -j^,Ps -1 + ^a - -^p. 



Appendix C: Physical quantities relevant for the 
estimate of the break frequency in blazars 
spectra 

C. 1. Redshift distribution of blazars 



In iMassardi et al.l (l2010al) the redshift distribution of radio 
sources at low frequencies is widely discussed, and we followed 
this paper to derive the redshift distribution of flat-spectrum 
sources. Most of the samples with redshift information do not 
distinguish between steep- and fl at-spectrum sources. Spectral 
information are present in the Kii hr et al.l(ll98lh catalogue, how- 
ever: this sample comprises 518 ERS to a 5-GHz flux density 
limit of 1 Jy, over an area of 9.81 1 sr. Based on the catalogued 
spectral indices, 299 sources are classified as flat-spectrum; 
212 of which are FSRQs (200 with measured redshift), 26 are 
BLLacs (20 with measured redshift) and 61 are classified as 
galaxies or with missing classification. Moreov er, in the Parke s 
quarter-Jy sample of flat-spectrum ERS (Jacks on"et al.l l2002h . 
reds hifts are available fo r the 58% of sources. From this sam- 
ple, iDe Zotti et al.l (|2005l) have defined a complete sub-sample 
of 514 objects with flux limit of 0.25 Jy, aiming at maximizing 
the fraction (~ 75%) of ERS with known redshift . This sub- 
sample includes 370 FSRQs (93% with known redshift) and 47 
BLLacs, of which only 21% with known redshift. 

In Fig. lC.ll we plot the redshift distributions of FSRQs and 
of BL Lacs, and the fits we used for our predictions. For FSRQs 
only, it has been possible to calculate them from both the sam- 
ples: in this case, the redshift distribution is observed to shift 
to higher redshifts as the flux limit of the sample is lowered 
down, with the peak of the distribution moving from z ^ 1.2 
to z ^ 1.5. Moreover, Fig. lC. II shows that the relative number 
of low-redshift FSRQs is strongly reduced, if a fainter flux de- 
tection limit is adopted. As for BLLacs, the redshift distribution 
can be obtained from the Kiihr et al. catalogue, and for only 20 
very bright objects. Because of the lack of information at faint 
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general, provides a better description of AGN jet properties, but 
has the disavantage to involve more free parameters than the ho- 
mogeneous model. Note, however, that in the case of the homo- 
geneous model it is assumed that all the observed X-ray flux is 
produced through inverse Compton scattering by the core com- 
ponent dominanting at the radio frequency. If part of the X-ray 
flux is produced in other components or by some other mech- 
anism, then the computed 5 is a lower limit. For these reasons 
and for greater simplicity, we assumed the same 6 distribution 
for BL Lacs and core-dominated quasars (in general for all the 
flat-spectrum sources), described by the fit in Fig. lC.ll 




Fig. C.l. Upper panels: redshift distributions of BLLacs (left 
panel) and of FSRQs (right panel) from the samples by 
iKuhr et al.1 (flQSl) (solid points) and Jackson et al. (2002) 
(empty points), and the corresponding fits (solid and dashed 
lines, respectively). Lower panels: h i stogra ms of the Doppler 
factor as obtained by iGhisellini et al.l (1 19931) (dashed lines) and 
IGu et alJ (1200 9^ (dotted lines) for BL Lacs (left panel) and core- 
dominated quasars (right panel). The solid line in the right panel 
is the fit of the two 6 distributions for core-dominated quasars. 



fluxes, this redshift distribution will be considered representative 
also for BL Lac sources with flux density lower than the sample 
limit. Owing to this, our predictions on number counts of ERS 
discussed in Sect. 7 are, therefore, more uncertain when applied 
to sources at 5 < 0. 1 Jy. 

C.2. The Doppler factor in AGNs 

The estimate of the Doppler factor 6 in AGNs is some- 
thing complex and model-dependent. I n the framewor k 
of the synchrotron s elf-Compton model (|Marscherl Il987|), 
iGhiseflini et al.l (Il993h calculated the Doppler factor for a sam- 
ple of 105 radio sources using VLBI measurements of the core 
angular dimensions and radio fluxes. The value of 6 is calcu- 
lated comparing the observed X-ray fluxes with the ones pre- 
dicted on the based of a homogeneous spherical emitting model. 
In Fig. lC.ll we show the 6 distribution for the 53 core-dominated 
quasars and the 33 BLLacs present in the sample. These dis- 
tributions can be compared with results from Gu et al. (2009) 
where the Doppler factor is computed using the Konigl inhomo- 
geneous jet model instead of the homogeneous spherical model. 
Their sample consists of 128 sources, with 80 core-dominated 
quasars and 26 BLLacs (37 quasars and 19 BLLacs are in com- 
mon with the sample used by Ghisellini et al.). The 5 distribu- 
tions are similar in the case of core-dominated quasars, with 
most of sources having 6 between 1 and 30, as expected for ob- 
jects where the relativisti c beamed emiss i on i s domin ant. For 
BLLacs the results from iGhiseUini et al.l (IT993) and Gu et al.l 
(12009 ) do not agree: the former find very low 6 values, extend- 
ing from 10"^ to 10; in the latter the 6 distribution is similar to 
the core-dominated quasars one. The inhomogeneous model, in 
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